arXiv: 1509.0228Ivl [astro-ph.EP] 8 Sep 2015 


Seasonal Variability of Saturn’s Tropospheric Temperatures, Winds and Para-H2 from 

Cassini Far-IR Spectroscopy 


Leigh N. Fletcher^*’'’, P.G.J. Irwin'’, R.K. Achterberg*^, G.S. Orton'*, F.M. Flasar*’ 

^Department of Physics Astronomy, University of Leicester, University Road, Leicester, LEI 7RH, UK 
^Atmospheric, Oceanic U Planetary Physics, Department of Physics, University of Oxford, Clarendon Laboratory, Parks Road, Oxford, 0X1 3PU, UK 
^Department of Astronomy, University of Maryland, College Park, MD 20742, USA. 

^Jet Propulsion Laboratory, California Institute of Technology, 4800 Oak Grove Drive, Pasadena, CA, 91109, USA 
^NASA/Goddard Space Flight Center, Greenbelt, Maryland, 20771, USA 


Abstract 

Far-IR 16-1000 spectra of Saturn’s hydrogen-helium continuum measured by Cassini’s Composite Infrared Spectrometer 
(CIRS) are inverted to construct a near-continuous record of upper tropospheric (70-700 mbar) temperatures and para-H2 fraction 
as a function of latitude, pressure and time for a third of a Saturnian year (2004-2014, from northern winter to northern spring). The 
thermal field reveals evidence of reversing summertime asymmetries superimposed onto the belt/zone structure. The temperature 
structure that is almost symmetric about the equator by 2014, with seasonal lag times that increase with depth and are qualitatively 
consistent with radiative climate models. Localised heating of the tropospheric hazes (100-250 mbar) create a distinct perturbation 
to the temperature profile that shifts in magnitude and location, declining in the autumn hemisphere and growing in the spring. 
Changes in the para-H2 {fp) distribution are subtle, with a 0.02-0.03 rise over the spring hemisphere (200-500 mbar) perturbed 
by (i) low-/p air advected by both the springtime storm of 2010 and equatorial upwelling; and (ii) subsidence of high-/p air at 
northern high latitudes, responsible for a developing north-south asymmetry in fp. Conversely, the shifting asymmetry in the para- 
H2 disequilibrium primarily reflects the changing temperature structure (and hence the equilibrium distribution of fp), rather than 
actual changes in fp induced by chemical conversion or transport. CIRS results interpolated to the same point in the seasonal 
cycle as re-analysed Voyager-1 observations (early northern spring) show qualitative consistency from year to year (i.e., the same 
tropospheric asymmetries in temperature and fp), with the exception of the tropical tropopause near the equatorial zones and belts, 
where downward propagation of a cool temperature anomaly associated with Saturn’s stratospheric oscillation could potentially 
perturb tropopause temperatures, para-H2 and winds. Quantitative differences between the Cassini and Voyager epochs suggest 
that the oscillation is not in phase with the seasonal cycle at these tropospheric depths (i.e., it should be described as quasi-periodic 
rather than ‘semi annual’). Variability in the zonal wind field derived from latitudinal thermal gradients is small (< 10 m/s per 
scale height near the tropopause) and mostly affects the broad retrograde jets, with the notable exception of large variability on the 
northern flank of the equatorial jet. The meridional potential vorticity (PV) gradient, and hence the ‘staircase of PV’ associated with 
spatial variations in the vigour of vertical mixing, has varied over the course of the mission but maintained its overall shape. PV 
gradients in latitude and altitude are used to estimate the atmospheric refractive index for the propagation of stationary planetary 
(Rossby) waves, predicting that such wave activity would be confined to regions of real refractivity (tropical regions plus bands at 
35 - 45'’ in both hemispheres). The penetration depth of these regions into the upper troposphere is temporally variable (potentially 
associated with stratification changes), whereas the latitudinal structure is largely unchanged over time (associated with the zonal 
jet system). 
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1. Introduction 

Far-infrared remote sensing (16-1000 jum) provides access to 
the three-dimensional structure of Saturn’s upper troposphere 
between approximately 70 and 700 mbar. This altitude range 
lies above the main condensible cloud layers (expected to be 
NH3 ice below the 1.2-bar level, Atreya et al., 1999), but 
contains both the radiative-convective (R-C) boundary and the 
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ubiquitous haze that is responsible for Saturn’s yellow-ochre 
hues. The R-C boundary signifies the transition from the 
convectively-unstable mid-troposphere to the statically-stable 
upper troposphere. The temperature field in this region is 
characterised by latitudinal contrasts between the cool zones 
(anticyclonic shear regions equatorward of prograde jets) and 
warmer belts (cyclonic shear regions poleward of prograde jets) 
(Conrath and Pirraglia, 1983). When Cassini first arrived at 
Saturn in 2004 (just after southern summer solstice at a plan- 
etocentric solar longitude of Ls - 293°), inversions of spec¬ 
tra from the Cassini Composite Infrared Spectrometer (CIRS, 
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Flasar et al., 2004) revealed stark hemispheric contrasts in up¬ 
per tropospheric temperatures and composition, caused by a 
combination of radiative, dynamical and chemical responses 
to the seasonal shift in sunlight (Fletcher et al., 2007). Qual¬ 
itatively similar asymmetries had been observed 24 years ear¬ 
lier by Voyager 1 (1980, - 8 . 6 °) and Voyager 2 (1981, 

Lj = 18.2°) (Flanel et al., 1981, 1982; Conrath and Pirraglia, 
1983; Conrath et al., 1998), although a direct quantitative com¬ 
parison of Cassini and Voyager findings only became possible 
in 2010-2011 when Saturn was at the same point in its seasonal 
cycle. 

Cassini has continued to study Saturn’s evolving atmosphere 
through spring equinox (L^ = 0°, August 2009) and on to¬ 
wards northern summer solstice (L^ = 90°, May 2017). Several 
studies have used CIRS 7-16 mid-infrared spectroscopy to 
monitor the evolution of asymmetries in temperature at pres¬ 
sures less than 250 mbar and stratospheric hydrocarbon com¬ 
position (Fletcher et al., 2010; Sinclair et al., 2013; Sylvestre 
et al., 2015; Fletcher et al., 2015). The broad temperature 
trends compare favourably to the predictions of Saturnian cli¬ 
mate models (Greathouse et al., 2010; Friedson and Moses, 
2012; Guerlet et al., 2014) that balance atmospheric heating 
from methane and aerosol absorption with cooling by hydro¬ 
carbon emission. This study utilises far-infrared (16-1000 /im) 
CIRS spectroscopy to investigate seasonal trends in the deeper 
troposphere, constructing a global record of tropospheric tem¬ 
perature and para-hydrogen variability over the ten-year span 
of the Cassini mission. 

Saturn’s 16-1000 //m spectrum is shaped by the collision- 
induced absorption of H 2 -H 2 and H 2 -He. Collisions between 
these molecules create instantaneous dipoles that form the con¬ 
tinuum, and the ratio of ortho-H 2 (corresponding to the odd ro¬ 
tational isomer of H 2 with parallel spins) to para-H 2 (the even 
isomer of H 2 with anti-parallel spins) can be deduced from the 
relative sizes of the broad S(0) (due to para-H 2 ) and S(l) (due 
to ortho-H 2 ) absorption features near 354 and 587 cm ', re¬ 
spectively. In Saturn’s deep troposphere the ratio between these 
two spin isomers is expected to be in equilibrium (3:1), equiva¬ 
lent to the para-hydrogen fraction {fp) of 0.25 (e.g., Massie and 
Hunten, 1982). CIRS constraints on fp mostly arise from the 
S(0) line and the translational continuum at wavenumbers less 
than 200 cm“', as the signal-to-noise ratio decreases substan¬ 
tially near the S(l) line. At the cold temperatures of the upper 
troposphere the equilibrium para-H 2 fraction (feqm) increases 
beyond 0.25 (reaching 0.4-0.5 at the tropopause) because of 
the increased population of para-H 2 energy levels at lower tem¬ 
peratures. However, if an upward-moving air parcel is rising 
quickly enough that the vertical mixing is faster than the chemi¬ 
cal equilibration timescale (i.e., assuming negligible ortho-para 
conversion during the motion), then the upwelling parcel will 
tend to retain its initial low para-H 2 fraction. The para-H 2 frac¬ 
tion would then appear to be in disequilibrium in Saturn’s upper 
troposphere. The spatial distribution of the para-H 2 fraction, 
and its deviation from equilibrium, can therefore be used to 
trace both vertical mixing and the efficiency of chemical equili¬ 
bration, with strong upwelling causing sub-equilibrium condi¬ 
tions { fp < feqm) and subsidence of cool tropopause air causing 


super-equilibrium conditions {fp > feqm) (Conrath et al., 1998; 
Fletcher et al., 2007). 

Voyager/IRIS revealed local sub-equilibrium conditions near 
60°S in early northern spring (Lj = 8.6 - 18.2° Conrath et al., 
1998), and an asymmetry with a higher fp in the spring hemi¬ 
sphere (super-equilibrium conditions from 70°N to 30° S). Con¬ 
versely, Cassini observations in late northern winter (Lj = 297- 
316°) revealed a region of low fp at the equator (interpreted 
as evidence of equatorial upwelling), super-equilibrium over 
much of the southern summer hemisphere and sub-equilibrium 
in the winter hemisphere poleward of 30°N (Fletcher et al., 
2007). The difference between the Voyager and Cassini results 
was believed to be seasonal, but a continuous Cassini time se¬ 
ries (rather than a single snapshot) is required to confirm this, 
and to distinguish between dynamical and chemical explana¬ 
tions for the observed fp distribution. 

The presence of aerosols may serve to catalyse the intercon¬ 
version between the ortho- and para-H 2 isomers (e.g., Massie 
and Hunten, 1982; Fouchet et al., 2003), so Saturn’s upper 
tropospheric haze is expected to be intimately linked to both 
the temperature and para-H 2 distribution. The composition of 
the haze is unknown, possibly linked to uplift of condensed 
NH 3 ice, subsidence of photochemically-produced hydrocar¬ 
bon smog, or photochemical destruction of tropospheric PH 3 
to form P 2 H 4 as a haze precursor (see the review by West et al., 
2009). The haze appears to be seasonally variable, and re¬ 
sponsible for the observed colour change from blue to yellow 
as Saturn’s northern hemisphere emerged from the depths of 
winter (e.g., Edgington et al., 2012). CIRS identified a local 
maximum in tropospheric temperatures (a ‘kink’ in the vertical 
temperature profile) between the tropopause and the radiative- 
convective boundary (Fletcher et al., 2007), that had been previ¬ 
ously identified at isolated southern latitudes by Voyager (Hanel 
et al., 1981; Lindal et al., 1985). The pressure level and mag¬ 
nitude of this temperature perturbation was latitudinally vari¬ 
able, being enhanced in southern summer and weak or absent 
in the northern winter hemisphere. The perturbation was also 
higher and weaker over the equator, and showed local maxima 
at ±15° latitude coincident with the neighbouring warm equato¬ 
rial belts. Fletcher et al. (2007) concluded that this was a radia¬ 
tive effect due to solar absorption Saturn’s upper tropospheric 
haze, and explained the asymmetry in terms of both seasonal in¬ 
solation and the variable distributions of tropospheric aerosols 
(later confirmed by radiative climate modelling by Friedson and 
Moses, 2012; Guerlet et al., 2014). This study attempts to track 
the evolution of the upper tropospheric temperature perturba¬ 
tion in relation to the distribution of hazes. 

Section 2 describes the processing and analysis of 
Cassini/CIRS far-infrared spectra to reconstruct a continuous 
record of tropospheric temperatures and para-H 2 fraction since 
2004 (Section 3). Section 4 then compares the results in early 
northern spring to those derived from Voyager/IRIS observa¬ 
tions to identify inter-annual variability, before discussing the 
implications of this study in Section 5. We will use the de¬ 
rived temperature and para-hydrogen fields to estimate the at¬ 
mospheric stratification and zonal wind field in three dimen¬ 
sions (altitude, latitude and time), calculating the implied po- 
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tential vorticity gradients and their implications for seasonally- 
variable mixing and the transmissivity of planetary waves in 
Saturn’s upper troposphere. In particular, measurements of the 
2D wind and potential vorticity gradients from Cassini data will 
provide insight into the dynamical processes occurring on the 
scale of Saturn’s jet system, a measure of the intrinsic stability 
of the features observed in the flow, and a closer connection be¬ 
tween observations and numerical simulations of giant planet 
atmospheres (Read et ah, 2006; Read et ah, 2009). 

2. Data Analysis 
2.7. Cassini/CIRS Data 

CIRS is made up of two interferometers sharing a com¬ 
mon telescope and scan mechanism (instrument and observa¬ 
tion strategies are reviewed by Flasar et al., 2004). This study 
primarily uses nadir spectra acquired with the longer wave¬ 
length (T > 16 jum) polarising interferometer (focal plane 1 ), 
measuring the 10-600 cm“^ region of Saturn’s far-infrared spec¬ 
trum using a single pixel with a circular held of view (4.1 
mrad diameter, with a Gaussian spatial response function of full 
width at half maximum of 2.5 mrad). Additional mid-infrared 
spectra sensitive to stratospheric methane emission (7.7 //m in 
focal plane 4) are included to constrain stratospheric tempera¬ 
ture. CIRS has a tuneable spectral resolution between 0.5 and 
15 cm ', and we use the 2.5 and 15 cm ' settings. The identifi¬ 
cation of trends in zonal mean atmospheric properties requires 
a significant number of points in a time series, so we follow 
Fletcher et al. (2015) and use both ‘targeted’ and ‘ride-along’ 
observations (i.e., those where CIRS acquires data during ob¬ 
servations by other Cassini instruments). The spatial and tem¬ 
poral distribution of the CIRS spectra used in this study are 
shown in Fig. 1. This incorporates both ‘sit-and-stare’ and 
scanning observations, and increases the temporal coverage be¬ 
yond the ‘FIRMAP’ datasets used in previous studies of Sat¬ 
urn’s temperature held (red crosses in Fig. 1 Fletcher et al., 
2007, 2010; Sinclair et al., 2013), at the expense of providing 
incomplete longitudinal coverage at a particular latitude and 
time. The spatial resolution of the observations were gener¬ 
ally greater for the 15-cm“^ data (2 - 5° latitude) than for the 
2.5 cm“* data (4-12°), although this was highly variable over 
Cassini’s orbital tour. 

All CIRS 2.5- and 15-cm ' nadir spectra from 2004 to 2014 
were extracted from the latest version of the CIRS database 
(version 4.3, Fig. 1) using month-long bins at monthly intervals 
and a 4° latitudinal grid (stepped every 2°), but applying strong 
Altering to omit spectra with calibration defects or inappropri¬ 
ate observing geometries, as described by Fletcher et al. (2015). 
Spectra were coadded within +10° of the mean zenith angle, 
and we required at least 40 target spectra and 4000 deep space 
calibration spectra for the coadded spectrum to be considered 
for spectral inversion. The coadded spectra were assumed to 
be representative of the zonal means, given that large-scale lon¬ 
gitudinal perturbations (e.g., thermal contrasts due to Saturn’s 
northern springtime storm, Fletcher et al., 2011; Achterberg 
et al., 2014) are extremely rare. Note that the storm-perturbed 
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Figure 1: Spatial and temporal distribution of more than 3000 CIRS focal plane 
one footprints at 2.5 cm“* spectral resolution (black dots) and 15.0 cm“* spec¬ 
tral resolution (red crosses). 

latitudes were omitted between 2011 and 2013 for the purpose 
of this study. This resulted in 1650 monthly-zonal-mean 2.5- 
cm“^ spectra and 1520 monthly-zonal-mean 15-cm ' spectra 
between 2004 and 2014. 

Examples of annual-mean far-IR spectra are shown in Fig. 
2 , where for the purposes of display the spectra were coadded 
over a full year within a 10 °-wide latitude bin and 10 °-wide 
zenith angle bin. Rotational lines of PH 3 and NH 3 are readily 
visible at 2.5 cm"* (Fig. 2c,d), but also cause oscillations of the 
50-200 cm"' spectrum at 15 cm"' resolution. Although formal 
spectral inversions are required to distinguish the contributions 
of temperature and the para-H 2 fraction, subtle changes in para- 
H 2 are manifest as changes to the spectral slope (e.g.. Fig. 2c), 
whereas shifts in temperatures generally scale the spectrum up 
and down (see Fig. 3). 

2.2. Spectral Inversion 

CIRS spectra for every month and latitude were inverted to 
derive atmospheric temperature and composition using the op¬ 
timal estimation retrieval algorithm, NEMESIS (Irwin et al., 
2008). A Levenburg-Marquardt iterative scheme is used to min¬ 
imise the residual between the forward model and the measure¬ 
ments, whilst maintaining physically-realistic solutions by us¬ 
ing smooth a priori state vectors (see Rodgers, 2000, for a full 
discussion of optimal solutions to the ill-posed inverse prob¬ 
lem). The a priori atmospheric state consists of: (a) a reference 
temperature profile based on previous CIRS analyses (Eletcher 
et al., 2007), which were themselves based on Voyager/IRIS 
and radio science measurements (Courtin et al., 1984; Lindal 
et al., 1985), defined on 120 pressure levels equally spaced in 
log p between 1 jubar and 10 bar; (b) the equilibrium para-H 2 
distribution; and (c) gaseous profiles for methane (Axed at a 
mole fraction of 4.7 x 10"^, Fletcher et al., 2009b), ammonia 
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Figure 2: Four examples of annually-averaged Cassini/CIRS far-lR (FPl) 
brightness temperature spectra, showing the difference between 15-cm“* res¬ 
olution (a, b) and 2.5-cm“* resolution (c, d). Spectra were averaged over a full 
year and a 10° latitude range, taking care to match emission angles. Not all lat¬ 
itudes are covered for every year. Warming trends in the northern hemisphere 
and cooling trends in the southern hemisphere are evident, as well as subtle 
shape changes due to altering para-H 2 . Error bars are plotted, but are smaller 
than the individual data points for wavenumbers below 450 cm“*. 


(Hurley et al., 2012), phosphine (based on the global mean re¬ 
sults of Fletcher et al., 2009a), and the stratospheric hydrocar¬ 
bons ethane and acetylene (using a spatial average of results 
from Guerlet et al., 2009). Rotational lines of NH 3 and PH 3 
have a non-negligible influence on the shape of Saturn’s spec¬ 
trum below 220 cm“^; stratospheric hydrocarbons form emis¬ 
sion features long ward of 650 cm-1, and so are included in 
our model for consistency only. 

Sources of spectral line data in our forward model are iden¬ 
tical to those in Table 4 of Fletcher et al. (2012). These absorp¬ 
tion coefficients are ranked according to their frequency distri¬ 
butions within a small spectral interval, producing smoothly- 
varying ^-distributions that increase the efficiency of transmis¬ 
sion calculations via the correlated-k approximation (Goody 
et al., 1989; Lads and Ginas, 1991). Absorption coefficients 
for the collision-induced H 2 -H 2 continuum were based on the 
results of Orton et al. (2007) (provided by M. Gustafsson for a 
range of fp values), who updated the calculations of Borysow 
et al. (1985). Additional opacities from H 2 -He, H 2 -CH 4 and 
CH 4 -CH 4 were included from Borysow and Frommhold (1987, 
1986); Borysow et al. (1988). These low-temperature opaci¬ 
ties are the same as those reported by Richard et al. (2012) 
for T < 400 K and included in the ‘Alternative’ folder of 
the HITRAN 2012 compilation (Rothman et al., 2013). The 
synthetic spectra were convolved with a triangular instrument 
function with a width appropriate to the spectral resolution of 
the measurements, approximating the Hamming apodisation of 
CIRS. Functional derivatives (or Jacobians, the rate of change 
of spectral radiance with the parameter of interest) calculated 
for Saturn’s reference temperature span 100-800 mbar in the 
10-500 cm“^ region and 1-3 mbar in the 1250-1350 cm“^ re¬ 
gion of stratospheric CH 4 emission (see Table 2 and Fig. 1 of 
Fletcher et al., 2007). Temperatures are smoothly interpolated 
between these two regions, and simply relax back to the a pri¬ 
ori for p > 800 mbar and p < 1 mbar. Para-H 2 sensitivity 
in CIRS far-IR spectra ranges from 100-550 mbar with a peak 
near 270 mbar, and CIRS cannot constrain fp at altitudes above 
the tropopause. 

For each spectrum, we simultaneously retrieved profiles of 
temperature and fp (in addition to scale factors for the a pri¬ 
ori distributions of PH 3 and NH 3 ) from the 50-550 cm * region 
of CIRS focal plane one and the 1230-1380 cm ' region of fo¬ 
cal plane four. The latter was added to the analysis to ensure 
that seasonal temperature contrasts in the lower stratosphere 
(Fletcher et al., 2010) were being correctly accounted for in 
the tropospheric retrieval. Fig. 3 shows why the temperature 
and para-H 2 are separable in the inversion, as they have dra¬ 
matically different effects on the shape of the collision-induced 
continuum. Temperature changes of ± 1 K modulate the contin¬ 
uum intensity in Fig. 3c, whereas changes in the para-H 2 frac¬ 
tion of +0.06 (the largest changes seen in this study) modify 
the spectral shape in Fig. 3b. These differences are easily de¬ 
tectable compared to the noise-equivalent spectral radiance on 
the coadded CIRS spectra. 

Although PH 3 and NH 3 affect the spectral shape between 50 
and 220 cm ', the retrieved mole fractions were found to be 
positively correlated with the retrieved para-H 2 fractions in the 
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upper troposphere because the fp retrievals use the translational 
continuum shortward of 200 cm ' in addition to the broad S(0) 
and S(l) lines (Fig. 3b). As a result, large excursions in one 
parameter would occasionally be mirrored by large excursions 
in the other. Rather than simply fixing the gases at their pri¬ 
ors, we allowed them to vary during the retrieval but filtered 
out substantial outliers from the final time series, as described 
in the next Section. Quantitative estimates of the para-H 2 frac¬ 
tion should be treated with caution, as they are sensitive to the 
helium mole fraction (e.g., Fletcher et al., 2007), which remains 
poorly constrained on Saturn. Our study fixes helium at a mole 
fraction of 0.135 as derived by Voyager/IRIS (Conrath and Gau¬ 
tier, 2000). CIRS focal plane 3 (600-1100 cm '), which also 
samples the upper troposphere, is omitted from this study be¬ 
cause the FPl and FP3 detector footprints don’t always overlap 
and therefore sample different horizontal regions. Although this 
is also true of focal plane 4 (used to constrain CH 4 emission), 
FP4 largely samples the stratosphere where spatial variability 
occurs over broader scales, and therefore has negligible effect 
on the tropospheric para-H 2 retrievals. Finally, it is interesting 
to note that a CIRS focal plane 2, which was descoped during a 
Cassini downsizing, would have provided enhanced signal-to- 
noise ratio near the S(l) line to improve our ability to constrain 
Saturn’s para-H 2 distribution. 

3. Cassini Results 

More than 3000 Tip) and fpip) profiles were derived from 
the CIRS spectra, representing zonal-mean conditions for each 
latitude and month. These are plotted in Fig. 4 to show the 
scatter of the measurements and the hemispheric asymmetries. 
Numerous outliers exist that could significantly corrupt the in¬ 
terpretation of a particular date and latitude - this could be due 
to (i) unexpected longitudinal variability so that the spectrum 
does not accurately represent a zonal average; (ii) calibration 
baseline offsets that are sometimes produced when the condi¬ 
tions vary between target CIRS spectra and the calibration spec¬ 
tra; or (iii) the combination of spectra observed from different 
orbital distances, and hence different fields of view averaging 
over spatial gradients. These spectra often appear to be real¬ 
istic, and are interpreted literally by the retrieval algorithm to 
produce a good fit. However, they produce retrieved tempera¬ 
tures and abundances that deviate noticeably from the general 
trends, and can be removed by filtering the results. We filtered 
out results where (i) Tip) and fpip) profiles exhibited substan¬ 
tial oscillations with altitude; (ii) the tropopause fp exceeded 
0.5; and (iii) the retrieved NH 3 and PH 3 scale factors deviated 
from the a priori by more than 2cr. This removed < 5% of the 
retrieved results. 

For the remaining data, the temporal coverage of a particular 
latitude was often sparse and incomplete (Fig. 1), so we fol¬ 
lowed Fletcher et al. (2015) in using smooth functions to inter¬ 
polate over the irregularly-spaced data to ‘reconstruct’ the tem¬ 
perature and para-H 2 variability. We required a quadratic func¬ 
tion to reproduce the 10 -year variability of temperatures at each 
altitude, whereas a linear function was found to be sufficient 
for the compositional variability. Examples of these temporal 
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Figure 3: Example of the effect on the CIRS spectral radiance of changing ei¬ 
ther the pai‘a-H 2 fraction or the tropospheric temperatures. The 2.5-cm“^ spec¬ 
trum shown as dots in panel a is for 45°N and 2014, as shown in brightness 
temperature units in Fig. 2c. The solid lines through these points represent 
a typical fit to the FPl spectmm. Panel a indicates that changes of order 10 
nW/cm^/sr/cm“* are barely discernible given the size of the radiances them¬ 
selves. However, the dilference plot in panel b shows that changing fp alters 
the shape of the spectral continuum, whereas panel c shows that changing the 
temperature alters the overall intensity of the spectrum. In b and c, the dotted 
line shows the CIRS noise-equivalent spectral radiance for a single spectrum; 
the dashed line shows the uncertainty when 100 spectra are coadded (a typical 
number for this study). 
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trends are shown for a selection of latitudes in Fig. 5. Note that 
these smooth functions are for interpolation only, and are not 
intended for extrapolation outside of the time period covered 
by the CIRS data. Clearly we have no knowledge of Saturn’s 
true temperatures during these temporal and spatial gaps, but 
we assume that the temperature and para-H 2 continue to follow 
the identihed trends. These gaps are fully accounted for in the 
uncertainty on the interpolation functions and the reconstructed 
fields. 

The interpolation functions were used to create ‘movies’ of 
the zonal-mean T{p) and fp{p) over the duration of the ob¬ 
servations, available in our Supplementary Online Material. 
Annually-averaged snapshots of the temperature are shown in 
Fig. 6; Fig. 8 shows temperature trends as a function of latitude 
and time at 100 mbar and 400 mbar; Fig. 9 shows the deriva¬ 
tive of the lapse rate {(fT Idz^) for temperature profiles shown 
in Fig. 10; and annual averages for the para-H 2 fraction and 
disequilibrium are shown in Fig. 11. 

Optimal estimation provides a formal estimate of the uncer¬ 
tainty on the retrieved profiles that depends on uncertainties in 
the measurements, a priori and forward model, and also on 
the degeneracies between the retrieved parameters. Examples 
of these individual uncertainties are shown in Fig. 5, and are 
used to weight the regression required to define the interpola¬ 
tion functions. Propagation of error from the retrieved quanti¬ 
ties to the derived quantities - the ‘reconstructed’ temperature 
and para-H 2 distributions in Section 3 and the winds, stratifica¬ 
tion, potential vorticity gradients and atmospheric refractivity 
in Section 5 - is a considerable challenge. Our approach is to 
use a Monte Carlo (MC) simulation to assess the potential range 
of these derived quantities. Each of the 3000-1- T{p) and fp{p) 
profiles was perturbed within the retrieval uncertainty envelope, 
using a Gaussian perturbation with (i) a peak altitude selected 
at random from a uniform distribution of pressures between 100 
and 800 mbar; and (ii) an amplitude calculated by multiplying 
the retrieval uncertainty by a number randomly selected from a 
normal distribution (i.e., the temperature profile remains within 
the Icr retrieval uncertainty range 68% of the time). We then re¬ 
calculated the regression lines for each time series to produce a 
new version of the reconstructed profiles, and repeat the deriva¬ 
tion of the thermal winds, buoyancy frequency, potential vor¬ 
ticity gradients and index of refraction (see Section 5). This 
process was repeated 100 times with random perturbations, and 
the resulting reconstructions are compared to provide a range 
of uncertainty in each parameter discussed in the following sec¬ 
tions. Note that this error analysis will not identify ‘unknown 
unknowns’, such as calibration biases or inadequacies in our 
forward model. 

3.1. Tropospheric Temperature Changes 

Saturn’s reconstructed tropospheric temperatures are shown 
in Eig. 6, , with the uncertainty estimated by the MC simulation 
shown in Eig. 7. Trends with time are shown in Eig. 8. The na¬ 
ture of uncertainties in regression lines means that the errors 
are time dependent, being smallest near the centre of the time 
series (near equinox in 2009) and larger towards the ends (2005 
and 2013). Uncertainties are larger at higher pressures where 
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Figure 4: Variability of tropospheric temperature and para-H 2 over the timespan 
of Cassini observations. The colour of each measurement represents the month 
of observation according to the key at the top of the figure. Thermal changes 
at 110 and 330 mbar (a,b) exceed the 0.5-1.0 K uncertainty on retrieved tem¬ 
peratures at each altitude. However, the scatter in retrieved fp from month to 
month and latitude to latitude (c,d) is considerably larger (and comparable to 
the uncertainty of 0.015) and masks para-H 2 trends at 300 mbar. 
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Figure 5: Examples of the polynomial fitting to temporal trends observed for temperature at 100 mbar (quadratic fits in panels a,b,c) and para-H 2 at 300 mbar (linear 
fits in panels d,e,f). In the top row, the comparison between 50°N and 50°S shows examples of seasonal warming and cooling, whereas the equator shows the effects 
of dynamics cooling the equatorial tropopause. The para-H 2 examples show a case where fp is increasing with time (d, 50°N); staying stable (f, 30°S); and an 
example where the linear fit is inappropriate (e, 30°N) near the location of the 2010-11 storm eruption (Achterberg et al., 2014). In this latter case, dashed lines 
show a two-step regression line to show how it deviates from the 10-year trend. 
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Figure 6: Snapshots of the zonal mean temperatures reconstructed from the 
sparsely-gridded temperature retrievals. Temperatures are averaged over a full 
year, and were chosen to represent Saturn’s equinox (2009) ±4 years. The full 
time-lapse sequence is available in the supplemental online material. 
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Figure 7: Temperature profile uncertainties in Kelvin calculated via an MC 
simulation in 2005, 2009 and 2013 for comparison with the changes observed 
in Fig. 6. Given the nature of uncertainties in a regression line, the eri'ors are 
time dependent, being smallest in the centre of the time series (2009). Contours 
are spaced every 0.1 K. 
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the sensitivity of the spectrum to temperature changes is lower, 
and at high latitudes where the gaps in temporal coverage are 
larger (Fig. 1). Indeed, uncertainties at northern high latitudes 
in 2005 are large (> 2 K) because these regions weren’t ob¬ 
served by CIRS until late 2006, and the interpolation functions 
are not meant to be used outside of these ranges. Furthermore, 
as Cassini typically observed from Saturn’s equatorial plane, 
the rings obscure the equator and result in gaps in the time se¬ 
ries (and larger uncertainties) at low latitudes. Over most of the 
planet the uncertainties are in the 0.4-0.8K range, as expected 
from the individual retrievals. 

Upon Cassini’s arrival at Saturn, tropopause temperatures 
showed a notable 10 + 2 K asymmetry between southern and 
northern high latitudes (Fig. 6a). This asymmetry has been re¬ 
versing throughout the duration of the mission, continuing the 
trend reported by Fletcher et al. (2010) using mid-infrared CIRS 
data. Southern hemisphere temperatures have generally cooled 
with the onset of autumn, although there are suggestions for a 
maximum in tropopause temperatures (92+2 K) in the 20-40°S 
region in 2006, a; 4 years after solstice (Fig. 8). In 2006, 100- 
mbar temperatures were less than 90 K only in the small region 
from 40-60°S, but by 2012 the entire southern hemisphere had 
cooled below 90 K, with temperatures less than 85 K poleward 
of45°Sby 2014. 

In the northern hemisphere. Fig. 8 shows the 100-mbar tem¬ 
peratures continuing to cool through 2007 and 2008, but then 
begin to warm as spring approaches. The lowest temperatures 
at the tropopause (78 + 1.5 K near 80 mbar) were detected sev¬ 
eral years after the winter solstice (2002). The precise timing 
of the temperature minimum varies as a function of altitude and 
latitude, with northernmost latitudes experiencing temperature 
minima near the spring equinox in 2009 (=« 7 years after sol¬ 
stice, consistent with Fletcher et al., 2015), whereas the min¬ 
imum occurred in 2007 at 45°N (=s 5 years after solstice. Fig. 
8 ). At higher pressures, the temperature minimum did not occur 
until early spring (a minimum north polar temperature of 87 + 1 
K at 220 mbar in 2010-2011 compared to 90 + 1 K by 2014), 
consistent with the increasing atmospheric lag times with depth. 
Indeed, Fig. 4b shows that 330-mbar temperatures poleward of 
75°N were relatively constant over the whole observation pe¬ 
riod. As of 2014 Fig. 4a shows that the northern tropopause 
was cooler than 85 K for all latitudes poleward of 45°N (with 
coolest tropopause temperatures of 82 ± 1 K in the north po¬ 
lar zone near 80-85°N), mirroring temperatures in the south¬ 
ern tropopause region poleward of 45°S. This suggests that the 
summertime asymmetry has now vanished and the two hemi¬ 
spheres show a distribution of tropopause temperatures that is 
approximately symmetric about the equator. 

The equatorial region, broadly defined as all latitudes 
equatorward of +20°, undergoes more unusual changes driven 
by dynamics. At 220 mbar, the 5 + 1 K contrast between the 
warm SEE and NEB observed at the start of the mission has 
now equalised (with 220-mbar temperatures of 95 + 0.5 K, 
in contrast to equatorial temperatures of 89 + 0.7 K, using 
uncertainties estimated by the MC technique). However, at 
higher altitudes the equatorial tropopause has cooled from 
89 + 0.5 K to 86 + 0.5 K over the ten year span of obser- 



2006 2008 2010 2012 2014 

Time 


Figure 8: Contour plot of temperature versus latitude and time for two tropo¬ 
spheric levels, 100 mbar and 400 mbar. These figures were generated from 
the polynomial-smoothed ‘reconstructed’ temperature fields, and are useful to 
show trends over the ten year span of the observations. Contours are spaced 
every 1.0 K for 100 mbar and every 0.75 K for 400 mbar to aid visualisation. 
Uncertainties on these temperatures are shown in Fig. 7. 


9 



















































vations (Fig. 8). The time-lapse movies suggest that this 
is due to the downward propagation of a cold temperature 
anomaly from the lower stratosphere, and such descent has 
been previously observed in the mid-stratosphere (0.01-20 
mbar) by authors studying Saturn’s semi-annual oscillation 
(SAO, e.g., Guerlet et al., 2011; Schinder et al., 201 1). Note 
that the air itself is not descending, but the pattern of warm 
and cool anomalies that result from complex wave-mean-flow 
momentum interactions that vary with time and space as the 
windshear changes. Indeed, the cool anomaly could be the 
result of local uplift and adiabatic expansion over a narrow 
altitude range. The variation in the equatorial tropopause 
suggests that this oscillating zonal-mean wind and temperature 
structure is also impacting the 50-150 mbar region. In the 
deeper atmosphere, the equatorial temperatures are stable with 
time. This zone remains one of the coolest regions compared 
to other latitudes, consistent with the idea of uplift and adia¬ 
batic expansion at low latitudes, mirrored by subsidence and 
warming over the neighbouring SEE and NEB. 

Einally, two anomalously warm northern hemisphere regions 
(T > 114K at 600 mbar) are visible in Eig. 6 in the deeper tro¬ 
posphere. The first region, the ‘bulge’ observed at 50 - 70°N in 
2005 that was identified in Cassini/CIRS inversions by Fletcher 
et al. (2007), cooled by 3-5 K over the ten year span of the data 
so that the latitudinal gradients in temperature {dT jdy) have 
vanished. The disappearance of this winter anomaly will have 
implications for changing potential vorticity gradients in Sec¬ 
tion 5. At the same time, a new region of T > 114 K developed 
in the 25 - 45‘’N region from 2010 onwards. This warming 
has been documented by Achterberg et al. (2014) and is asso¬ 
ciated with the eruption of Saturn’s northern springtime storm. 
However, a similar warm region was noted one Saturn year ear¬ 
lier by Hanel et al. (1982) and Conrath and Pirraglia (1983), 
and we speculate that this deep atmospheric warming is a dy¬ 
namical response to the changing season. There appear to be 
no comparable regions in the southern hemisphere during this 
time period, with 600-mbar contours remaining relatively uni¬ 
form with latitude. 

3.2. Aerosol Heating Effect 

The region of localised heating in the 100-250 mbar region 
(Section 1) can be quantitatively investigated via measurements 
of the ‘curvature’ of the T(p) profile by taking the second 
derivative, dTjdz - d^T jdz^ (where F is the lapse rate). This 
thermal perturbation is believed to be related to localised heat¬ 
ing within Saturn’s upper tropospheric haze (Fletcher et al., 
2007; Friedson and Moses, 2012; Guerlet et al., 2014), and 
showed a pronounced asymmetry early in the mission with a 
stronger perturbation in the southern summer hemisphere, and 
negligible in the north. The magnitude of the curvature is shown 
in Fig. 9, and the reconstructed T{p) profiles for three selected 
latitudes are provided in Fig. 10. Uncertainties tend to be mag¬ 
nified wherever derivatives are involved, so the MC simulation 
is used to show the uncertainty on this curvature in the bot¬ 
tom row of Fig. 9, in units of a single contour spacing (0.002 
K/km^). The uncertainty in the curvature is typically less than 
10% of a contour interval, but can rise to 25% and higher at the 


equator, poles and the storm perturbed latitudes due to gaps in 
the time series. The peak negative curvature (grey shading) is 
located at the pressure of maximum thermal perturbation, and 
the peak positive curvature is located where the T{p) profile 
tends to the deeper adiabat. The inflection points (where the 
temperature profile changes the direction of curvature) are used 
to provide approximate upper and lower pressure limits on the 
thermal perturbation. Note that regions of negative curvature 
near 600 mbar are likely to be spurious, related to relaxation of 
our retrieved T{p) back to the a priori profile. Finally, the size 
of the curvature is sensitive to our choice of vertical smooth¬ 
ing in the original inversions (which is kept the same for every 
retrieval, see Fletcher et al., 2007, for details), so we restrict 
our discussion below to qualitative trends in the curvature over 
time. 

During southern summer, the largest thermal perturbation is 
observed over the SEE, and the feature is observed over much 
of the southern hemisphere and the NEB in the 150-250 mbar 
range. The perturbation appears higher (70-150 mbar) and 
weaker over the equator, as observed by Eletcher et al. (2007), 
although uncertainties are highest at low latitudes. There were 
no regions of negative curvature in the northern winter hemi¬ 
sphere poleward of 20°N. The magnitude of the perturbation 
shrank but the location remained unchanged through to the 
spring equinox in 2009. Erom 2008 onwards, regions of nega¬ 
tive curvature began to appear in the northern hemisphere, start¬ 
ing in the cool zone near 45°N. By 2010-11, the thermal pertur¬ 
bation was being detected throughout the northern hemisphere 
poleward of 35 - 40°N. At the same time, the curvature had 
shrunk so much in the south that it vanished entirely from lat¬ 
itudes poleward of 45 - 50‘’S, giving the impression that the 
asymmetry observed at the start of Cassini’s mission had re¬ 
versed. By 2014, the NEB and SEB exhibited thermal perturba¬ 
tions of the same relative size, and the ‘kink’ was more preva¬ 
lent in the north than in the south - a change larger than the 
uncertainties shown in Eig. 9. Interestingly, the high extension 
over the equatorial zone vanished over the ten-year span of the 
data (the peak negative curvature was pushed deeper from near 
80 mbar to near 150 mbar), which could be related to the cool¬ 
ing equatorial tropopause in response to Saturn’s SAO wave 
train. 

With CIRS data alone, it is impossible to assess whether the 
growth of the thermal perturbation in the north is due simply 
to enhanced sunlight on an existing, static tropospheric haze 
(as suggested by the circulation and transport model of Eried- 
son and Moses, 2012), or whether the dominant effect is due to 
changes in the haze itself as spring progresses (e.g. Karkoschka 
and Tomasko, 2005). The latter is supported by colour changes 
observed by Cassini’s cameras (Edgington et al., 2012), and a 
stark asymmetry in tropospheric haze opacity in 2006 was iden¬ 
tified by Cassini’s near-infrared spectrometer (VIMS, Baines 
et al., 2006), with negligible haze opacity in the hemisphere 
emerging from winter. Growth of tropospheric aerosols in 
springtime was suggested by Karkoschka and Tomasko (2005) 
from their analysis of Hubble data, and Guerlet et al. (2014) 
show that the size of the thermal perturbation would increase 
with tropospheric particle size and/or increased opacity. How- 
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Figure 10: Examples of the reconstructed T{p) profiles showing the location of 
the ‘kink’ related to aerosol heating in Saturn’s upper troposphere. One curve 
is shown for each of the ten years. Thi'ee latitudes (the 60°N, the equator and 
15°S) are shown, with 10 and 20-K offsets applied to the latter two. Blue curves 
indicate profiles from early in the mission, red curves indicate curves from later 
in the mission, showing how the curvature of the T{p) profiles has altered over 
time. 

ever, a quantitative time series of Saturn’s haze variability is 
required to fully understand the shifting location and size of 
this thermal perturbation. 

3.3. Para-Hydro gen 

Zonal-mean cross-sections for para-H 2 as a function of time 
were generated from the linear interpolation function, and five 
annual snapshots of the time-lapse movies are shown in Fig. 
11. As we found for the temperature field, uncertainties on fp 
estimated by the MC simulation (Fig. 12) are time-dependent, 
being smallest near the middle of the time series (=s 0.005 for 
300 mbar in 2009) and growing at either end. Fig. 5 indicated 
that a single linear rate was not appropriate at the storm lat¬ 
itudes in 2011-2013 (15 - 45°N), where fp exhibited a sharp 
drop (Achterberg et al., 2014) (see Fig. 1 1(d)), so we use a 2- 
step linear function (pre-2011 and post-2011) for northern mid¬ 
latitudes. These linear rates of change of para-H 2 {d{fp)ldt) 
and their uncertainties are shown in Fig. 13, broken down into 
three altitude ranges. The smaller number of points contribut¬ 
ing to the regression lines post-2011 result in larger uncertain¬ 
ties (> 0.015) in the reconstructed fp at the storm latitudes in 
Fig. 12c. Furthermore, we caution that CIRS is only sensitive to 
para-H 2 in the 100-550 mbar range (Fletcher et al., 2007), and 
cannot constrain f, at altitudes above the tropopause. Finally, 
we remind the reader that the quantitative values of fp and feqm 
are sensitive to Saturn’s helium abundance, for which we use 
the value of Conrath and Gautier (2000). For these reasons we 
restrict our discussion to changes to the para-H 2 fraction with 
time, rather than the absolute values. 

Fig. 13 shows that the majority of the retrievals at 100 mbar 
are consistent with no change to fp over the ten-year data span 
for much of the mission (uncertainties on d(fp)ldt mean that 


Planetographic Latitude 

-90-75-60-45 -30 -15 0 15 30 45 60 75 90 
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Figure 12: Uncertainties on the reconstructed para-hydrogen fraction calculated 
from an MC simulation within the formal retrieval envelope. The uncertainty 
is time dependent, and annually-averaged examples are shown for 2005, 2009 
and 2013. Contours are spaced every 0.001, with grey shading used to indicate 
where uncertainties exceed 0.015 (equivalent to the contour spacing in Fig. 11), 
particularly following the 2010-2011 storm. Uncertainties on the fp disequilib¬ 
rium show the same spatial structure and similar amplitudes. 
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Figure 9: The curvature of the temperature profile as a function of time, used as a proxy for the size of the aerosol-related temperature perturbation in Saturn’s 
upper troposphere. Contour intervals are spaced by 0.002 K/km^, with grey shading used to indicate regions of strong negative curvature (i.e., the strongest thermal 
perturbation). Dotted contours indicate regions of positive curvature. Contours are shown as annual means for the equinox ±5 years, with the full time-lapse 
sequence available in the supplementary materials. Regions of negative curvature at p > 500 mbar are spurious, caused by the retrieved Tip) relaxing to the prior in 
a region of low spectral information content. The lower panel shows the estimated uncertainty on the curvature, shown as a percentage of the 0.002 K/km^ contour 
spacing (i.e., typical values of 10-25% of a contour interval. 


they are consistent with zero). There is a hint that fp has in¬ 
creased over the southern hemisphere by 0.01-0.02, but this is 
at the level of the uncertainty described above. The only signif¬ 
icant change at 100 mbar is in the fp poleward of 45°N, where 
para-H 2 has increased by 0.03 ± 0.01 over the course of the 
mission. At intermediate pressures (200 mbar) near the peak of 
the contribution function, the southern hemisphere increase of 
0.01 is still apparent, as is the 0.02-0.03 increase poleward of 
45°N. Before 2011, fp had shown an increase with time near 
30°N (Fig. 5e and the red crosses in Fig. 13), but fp decreased 
abruptly between 15 - 45°N at the start of 2011 in associa¬ 
tion with upwelling during the springtime storm (Achterberg 
et ah, 2014). In the aftermath of the storm, fp continued to rise 
but at a faster rate than before (blue squares in Fig. 13). At 
higher pressures (500 mbar) there is no change over the south¬ 
ern hemisphere, a 0.01-0.02 increase in fp at the equator, a 
general 0.02-0.03 increase over northern latitudes (particularly 
50 - 70°N) punctuated by effects of the northern storm, and 
negligible change over the north polar region. 

Fig. 11 and 13 reveal that Saturn’s tropospheric para-H 2 frac¬ 
tion has generally increased over the northern hemisphere, but 
that this trend is modified in both the storm region and the 
northern high-latitudes. The latter can be interpreted as high- 
latitude tropospheric subsidence in the spring season, consis¬ 
tent with similar subsidence identified in the stratosphere by 
Fletcher et al. (2015). At higher altitudes, Fletcher et al. (2015) 
used the rising north polar temperatures and enhancements in 


hydrocarbons ethane and acetylene to infer downward motion 
over the 75-90°N region. This subsidence may persist into the 
upper troposphere for pressures less than 300-400 mbar, but not 
as deep as 500 mbar where fp was not observed to be chang¬ 
ing. Furthermore, analysis of the Voyager/IRIS data by Con- 
rath et al. (1998) also suggested a higher fp in the northern 
hemisphere just after spring equinox, and suggested that this 
was due to a downward displacement. They found large values 
of fp centred on 60°N (confirmed by our reanalysis in Section 
4), exactly where Fig. 13c shows some of the largest temporal 
trends at 500 mbar outside of the storm latitudes. The increase 
we see at 60°N has been continuous over Cassini’s observa¬ 
tions, unperturbed by the effects of the storm. The changes in 
this region, and in the north polar region, are likely to be due 
to responses of the fp distribution to dynamics. Conversely, 
the small 0.01-0.02 rise in fp at the tropopause over the south¬ 
ern hemisphere as autumn approaches could be a seasonal re¬ 
sponse related to shorter equilibration times in the hazier south¬ 
ern hemisphere. 

The changing para-H 2 disequilibrium observed in Fig. 11 re¬ 
flects primarily the changing temperature structure (and there¬ 
fore distribution of fp,eqm) rather than actual changes in fp 
induced by chemical conversion or transport. Uncertainties 
on the disequilibrium are slightly larger than those shown in 
Fig. 12 for fp alone, as they must account for both temper¬ 
ature and fp uncertainties. Nevertheless, they are generally 
smaller than 0.01 (the contour spacing of the disequilibrium 


12 


Pressure [bar] Pressure [bar] 




















































Planetographic Latitude Planetographic Latitude 



-90 -75 -60 -45 -30 -15 0 15 30 45 60 75 90 -75 -60 -45 -30 -15 0 15 30 45 60 75 90 
Planetographic Latitude Planetographic Latitude 


Figure 11: Reconstructed para-H 2 fraction (left column) and degree of disequilibrium (right column) for five years: Saturn’s northern spring equinox (2009) and 
four years either side of it. Uncertainties calculated via our MC simulation are shown in Fig. 12. Contours of fp are spaced every 0.015, equivalent to the largest 
uncertainty. Contours of disequilibrium are spaced every 0.01, with sub-equilibrium conditions {fp < feqm) indicated by the grey shading and dotted contours. The 
full time-lapse sequence is available in the supplemental online material. 
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Figure 13: Rates of change of para-H 2 fraction at three different altitudes: (a) 
100 mbar; (b) 200 mbar; (c) 500 mbar. Black circles use a linear regression 
over the whole span of the data, red crosses between 15 - 45°N use only pre- 
2011 data for the regression to eliminate the springtime storm effects; blue 
squai'es use only post-2011 data to show the increased rates in the aftermath of 
the storm. The rising fp at high northern latitudes is more apparent at lower 
pressures. The horizontal dotted line shows where fp has been invariant with 
time. 


plots in Fig. 11. At the start of the mission, Fletcher et al. 
(2007) identihed super-equilbrium conditions { fp > fp^eqm) over 
the southern summer hemisphere and sub-equilibrium condi¬ 
tions {fp < fp^eqm at the equator and over the northern win¬ 
ter hemisphere. They speculated that this difference was not 
due to large-scale inter-hemispheric circulation (i.e., upwelling 
in the north and subsidence in the south), but rather due to 
the enhanced efficiency of conversion of para-H 2 to ortho-H 2 
in the summer hemisphere where there was a greater density 
of aerosols to serve as catalysts (e.g., Fouchet et al., 2003). 
Studying the snapshots of the disequilibrium in Fig. 11, we 
see that neither of these interpretations is correct. The super¬ 
equilibrium conditions have disappeared in the south as the 
temperatures dropped, tending towards equilibrium and sub¬ 
equilibrium (poleward of 45°S in 2014) as fp^eqm rose to match 
the value of fp that was detected in southern summer. Likewise, 
in the north the sub-equilibrium conditions have declined as 
the temperature rose, so that equilibrium and super-equilibrium 
conditions are detectable between 45 - 75°N for p > 200 mbar. 
In this case, fp^eqm decreased with the warming springtime tem¬ 
peratures to match the fp found there during northern winter. 
In short, the large-scale asymmetries in the disequilibrium are 
simply driven by the temperature-dependent changes to fp^eqmi 
and not due to changes in fp itself. Neither large-scale circula¬ 
tion nor aerosol catalysis are required to explain the asymmetry 

fp ~ fp,eqm- 

At low pressures, the detection of disequilibrium at certain 
points in Saturn’s seasonal cycle is consistent with the find¬ 
ings of Conrath et al. (1998), whereby the para-H 2 equilibration 
timescale greatly exceeds the timescales for thermal changes 
(radiative or dynamical), so the two spin isomers are ‘frozen in’ 
and fp (but not fp - fp^eqm) traces the atmospheric circulation. 
However, fp generally remains close to equilibrium at higher 
pressures, and we cannot rule out the possibility that surface 
catalysis on Saturn’s aerosols plays a role in rapid interconver¬ 
sion of the two spin isomers (e.g., Massie and Hunten, 1982). 
As the timescales for equilibration (frozen or fast) influence the 
lapse rate, this degeneracy will influence the quantities derived 
in Section 5. 

The pattern of shifting disequilibrium is only broken where 
regional dynamics become important at the equator, poles and 
storm-perturbed latitudes. The sub-equilibrium conditions at 
the equator and corresponding super-equilibrium conditions 
near +17° have persisted throughout the dataset, although there 
are signs that the local minimum in fp has weakened at higher 
pressures (e.g.. Fig. 13c). This pattern was interpreted at the 
start of Cassini’s mission as due to upwelling in the equato¬ 
rial zone and subsidence in the neighbouring equatorial belts 
(e.g., Fletcher et al., 2007). The balance between upwelling 
and chemical equilibration could now be shifting, possibly 
due to weakened (or reversing) vertical motions influenced by 
the oscillating temperatures and winds at lower pressures. At 
other locations the storm eruption caused fp to fall between 
15 - 45°N to generate a corresponding sub-equilibrium region 
in Fig. 1 1(d). And at high northern latitudes, where the largest 
sub-equilibrium was revealed a the start of the Cassini mission, 
the high-latitude subsidence is causing the most rapid return to 
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equilibrium conditions anywhere on the planet (a change in fp 
of 0.04-0.05). 

4. Voyager/IRIS Reanalysis 

The reconstructed Cassini T{p) and fpip) span the same 
points in the seasonal cycle (L^ = 8.6°, April 2010, and 
Lj = 18.2°, February 2011) as the Voyager 1 and 2 flybys, re¬ 
spectively, and therefore permit a quantitative comparison of 
tropospheric conditions to search for evidence of inter-annual 
variability. This has previously been investigated by Li et al. 
(2013) and Sinclair et al. (2014), but using Cassini and Voyager 
datasets which were not obtained at exactly the same point in 
Saturn’s season. The Voyagers carried the Infrared Radiome¬ 
ter, Interferometer and Spectrometer (IRIS) instruments, with 
a Michelson interferometer providing 180-2500 cm ' spectra 
of the giant planets at a spectral resolution of 4.3 cm ' (Hanel 
et al., 1980). We reanalyse the Voyager/IRIS Saturn data (Con- 
rath et al., 1998) to ensure that the same spectral model (in 
particular, the updated H 2 -H 2 collision induced absorption), a 
priori atmosphere and retrieval assumptions are used for both 
Voyager and Cassini datasets. IRIS spectra were extracted from 
the expanded volumes available on NASA’s Planetary Data Sys¬ 
tem', and zonal-mean spectra were constructed from observa¬ 
tions within +2 days of closest approach (November 13, 1980 
and August 26, 1981 for Voyager 1 and 2, respectively). Unlike 
Cassini’s global mapping sequences, the Voyager data were ac¬ 
quired in individual latitude scans resulting in irregular latitu¬ 
dinal coverage and low spatial resolutions at mid-to-high lati¬ 
tudes. Zonal-mean spectra were coadded in latitudinal bins 8° 
wide and stepped every 2°, ensuring that the IRIS field of view 
was entirely on Saturn’s disc. Spectra were pre-filtered to re¬ 
move any with a high degree of oscillations, and uncertainties 
were derived from the Voyager noise-equivalent spectral radi¬ 
ance (NESR), as described by Sinclair et al. (2014). In total, 
we used 990 Voyager-1 spectra (covering all latitudes) out of 
approximately 2100 spectra taken within ±2 days of closest ap¬ 
proach. Voyager-2 spectra had a limited latitudinal coverage 
(covering the northern hemisphere, with limited observations 
in the south) and proved more difficult to fit, possibly due to 
a misalignment of the interferometer (Hanel et al., 1982; Sin¬ 
clair et al., 2014). Furthermore, their Cassini-era counterpart 
coincides with a time where Saturn’s northern latitudes were 
strongly perturbed by the 2010-2011 storm, so they were omit¬ 
ted from the remainder of this study. 

We derive T{p) and fpip) from the 220-550 cm“' region 
of the IRIS spectra (encompassing the S(0) line and the long- 
wavelength edge of S(l)), supplemented by the 1230-1370 
cm ' region to provide constraints on stratospheric tempera¬ 
tures (Sinclair et al., 2014). The improved signal-to-noise ratio 
of IRIS near the S(l) line compared to CIRS, but the absence 
of data in the translational part of the spectrum below 220 cm ' 
(e.g.. Fig. 3), implies that IRIS offers different spectral con¬ 
straints to CIRS for para-H 2 , so caution must be taken in quan¬ 
titative comparisons. Even if both instruments were restricted 


'http://pds-rings.seti.org/voyager/iris/expanded_volumes.html 


to an identical spectral range, the different noise characteristics 
of IRIS versus CIRS would inevitably provide different mea¬ 
surement uncertainties in the retrieval process. 

Erom inversions of the Voyager-1 spectra, we confirm that 
Saturn’s 1-mbar equatorial temperatures were colder than the 
neighbouring belts at this time period, supporting the conclu¬ 
sion of Sinclair et al. (2014) that Voyager and Cassini captured 
a different phase of Saturn’s stratospheric oscillation (known as 
the semi-annual oscillation, SAO). Tropospheric temperatures, 
the thermal perturbation due to aerosol heating, para-H 2 and the 
degree of disequilibrium are shown side by side for Voyager-1 
and Cassini in Eig. 14, with the Cassini results interpolated to 
April 2010 to match Lj = 8.6°. Typical errors on the IRIS T(p) 
are 0.6-0.7 K between 100 and 600 mbar, 1.6-1.8 K at 1 mbar. 
Typical errors on the fp(p) are 0.015-0.019 between 100-500 
mbar. These are comparable to the uncertainties on the CIRS 
reconstructions in Eig. 6 and Eig. 11 . 

Tropopause temperatures in Eig. 14a show a similar asym¬ 
metry for both datasets between the cool north and warmer 
southern hemisphere, although IRIS found cooler (85-87 
K) southern hemisphere temperatures than CIRS (88-90 K), 
whereas the northern tropopause temperatures (80-85 K) were 
roughly comparable between both epochs. At pressures of 500 
mbar and deeper, the IRIS inversions also produce systemati¬ 
cally cooler temperatures by 3-4 K compared to the CIRS in¬ 
versions. Random uncertainties from the IRIS inversion and 
CIRS interpolation procedure are smaller than 1 K at these alti¬ 
tudes, so cannot explain this offset. This could be a real change 
over the Saturnian year, but systematic calibration offsets be¬ 
tween the two instruments have not been adequately assessed, 
and these would significantly hamper any attempt at quantita¬ 
tive comparison. We therefore restrict our comparison to rela¬ 
tive spatial contrasts. 

The largest differences between the two epochs occur within 
15° the equator; in the Voyager-1 observations the NEB ap¬ 
pears warmer than the SEB by 2-3 K at 100 mbar, whereas the 
Cassini observations show the NEB to be cooler than the SEB 
by 1-2 K, qualitatively consistent with the results from Li et al. 
(2013). This difference is only apparent for p < 200 mbar, but 
serves to alter the latitudinal thermal gradient around the equa¬ 
torial zone and therefore the vertical shear on the zonal wind 
(Li et al., 2013). The temperature comparison therefore sug¬ 
gests non-seasonal effects at Saturn’s equatorial tropopause that 
disrupts the seasonal cycle. Also notable is the 20 - 45°N warm 
region in the IRIS retrieval at 600 mbar (also evident in Eig. 6 
of Conrath et al., 1998) that is absent from the CIRS retrieval 
(this was 8 months prior to the eruption of the springtime storm 
in December 2010). 

Cassini’s results showed that the thermal perturbation be¬ 
tween 150-250 mbar featured a maximum curvature elevated 
over the equator, with the size of the perturbation decreasing in 
the southern autumn hemisphere and growing in the northern 
spring hemisphere. However, the perturbation appeared larger 
in the north during the Voyager-1/IRIS observations than dur¬ 
ing the Cassini/CIRS observations, suggesting that the mech¬ 
anism responsible for the heating may not be exactly repro¬ 
ducible from year to year. Unfortunately, detailed comparisons 
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(a) Temperature [K] 
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Figure 14: Comparison of results from Voyager/IRIS (November 1980) with those from Cassini/CIRS (April 2010), interpolated to the same point in Saturn’s 
seasonal cycle (L^ = 8.6°). Panel (a) shows tropospheric temperatures with a 2-K contour spacing; panel (b) shows the para-H 2 fraction with a 0.015 contour 
spacing; panel (c) shows the deviation of fp from equilibrium in steps of 0.01; and panel (d) shows the curvature of the temperature profile in steps of 0.002 K/km^. 
Despite the similar season, stark differences are evident in the tropical regions and at mid northern latitudes, although these are close to the conservative uncertainty 
estimates (1 K on T{p) and 0.02 on fp(p)) for the two datasets. 
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of the aerosol loading between the Voyager and Cassini spring¬ 
time epochs are not yet available to see if this ties to the thick¬ 
ness and coloration of Saturn’s tropospheric hazes. However, 
we caution that quantitative comparisons of the T{p) curvature 
between the two datasets are difficult, as the lower-resolution 
IRIS data provides a different constraint on the spectral inver¬ 
sions compared to the higher-resolution CIRS data. 

Finally, the para-Ha fraction shows an asymmetry between 
the northern and southern hemispheres in the Voyager epoch, 
whereas the Cassini results are more symmetric with a gen¬ 
erally higher value of fp in the northern spring hemisphere 
(by 0.02-0.03). There are, however, several quantitative dif¬ 
ferences between the two epochs: Voyager observed the high¬ 
est fp near 60°N that is not apparent in the CIRS inversions; 
the hemispheric asymmetry appears more pronounced in the 
IRIS inversion; and the equatorial minimum in fp appears more 
subdued in the IRIS inversion than in the CIRS results. Note 
that these quantitative differences are within the « 0.02 un¬ 
certainty on the CIRS and IRIS inversions, and could be ne¬ 
glecting important systematic calibration differences between 
the two datasets. The degree of disequilibrium in Fig. 14c 
shows this asymmetry more vividly - sub-equilibrium condi¬ 
tions cover much of the northern hemisphere in the CIRS epoch, 
whereas super-equilibrium conditions are apparent at p > 200 
mbar from 15 - 60°N in the Voyager epoch. Furthermore, 
the sub-equilibrium conditions in the equatorial zone do not 
extend as high during the Voyager observations as during the 
Cassini observations. However, the contours in Fig. 14c are 
spaced every 0.01, similar to the uncertainty on fp, so these 
apparent offsets are may not be particularly significant. Dif¬ 
ferences at the equatorial tropopause may be related to inter¬ 
annual variability in the oscillating wind and temperature fields 
above the tropopause, whereas the different hemispheric con¬ 
trasts are harder to explain. The differing spectral constraints 
offered by CIRS and IRIS may also be corrupting the compar¬ 
ison, but we cannot entirely discount a difference in the phas¬ 
ing of the seasonal inter-hemispheric transport between the two 
epochs (i.e., the strength of the sinking motion at springtime 
high latitudes). 

In summary, although hemispheric contrasts in temperature 
appear to be broadly reproducible from one Saturn year to the 
next, hemispheric asymmetries in the para-H 2 fraction appear 
different between the two datasets. Voyager-1 found the largest 
fp at northern high latitudes, and it is notable that Cassini has 
identified an increasing fp at the same location as springtime 
progressed. Equatorial conditions are modulated by a dynamic 
perturbation that is not locked to the seasonal cycle. As sug¬ 
gested by the stratospheric results of Sinclair et al. (2014), we 
propose that Saturn’s equatorial oscillation has a period that dif¬ 
fers from half a Saturn year at the tropopause (i.e., it is not 
rigorously ‘semi-annual’ and should be referred to as a quasi- 
periodic oscillation). 

5. Discussion 

Cassini/CIRS provides the capability to measure the zonal 
mean temperature structure as a function of time, latitude and 


pressure, from which we can derive several atmospheric diag¬ 
nostics to understand the dynamic response to seasonal vari¬ 
ability. These include: (i) the static stability of the atmosphere 
as a measure of its susceptibility to convective motions; (ii) use 
of the thermal windshear relation to determine the zonal wind 
field as a function of altitude, and therefore extending studies 
of vorticity gradients far above the cloud-tops; and (iii) the im¬ 
portance of these vorticity gradients in confining and guiding 
planetary wave activity (e.g. Achterberg and Flasar, 1996; Li 
et al., 2006). Given a suitable and accurate radiative-convective 
model, deviations of the thermal structure from the expecta¬ 
tions of radiative balance could be used to estimate the sense 
and magnitude of vertical motions (e.g. Conrath et al., 1990), 
but we leave that to a future study given the difficulties in repro¬ 
ducing Saturn’s tropospheric temperature and aerosol response 
as a function of season (Guerlet et al., 2014). We explore some 
of these quantities, derived directly from the thermal field, in 
the following sections. Uncertainties in these derived fields use 
the MC simulation described in Section 3. 


5.7. Atmospheric Stratification 

The seasonal change in the upper troposphere might be ex¬ 
pected to alter the stability of the atmosphere to convective 
motions, by changing the relative altitude of the radiative- 
convective (R-C) boundary over time. This can be studied via 
the frequency of buoyancy oscillations {N, known as the Brunt 
Vdisala frequency): 




T\dz Cpj 


( 1 ) 


Here g is the gravitational acceleration calculated as a function 
of altitude and latitude accounting for Saturn’s oblateness; Cp 
is the specific heat capacity calculated for the specific tempera¬ 
ture and para-H 2 fraction at each altitude and latitude; the quan¬ 
tity Fa - g/cp is the dry adiabatic lapse rate and F = -dTjdz 
is the lapse rate measured from the reconstructed temperature 
profiles. N allows us to compare the stratification from year to 
year directly - statically unstable regions have < 0 so that 

N is imaginary; statically stable regions have > 0 so that 

N is real and represents a simple harmonic motion of buoyancy 
oscillations. 

The uncertainty on N can arise from numerous sources - the 
reconstructed temperature field appears as both a divisor and 
a derivative in the equation above, and both the temperature 
and para-H 2 fraction determine the specific heat capacity as a 
function of altitude and latitude. The MC simulation shows the 
range of values of N expected from the propagation of uncer¬ 
tainties in T and fp. The error range on the 2009 values of 
N (Fig. 15) varies from 1 x 10^^ s“* at mid-latitudes and 200 
mbar to 1.2 x 10~^ s“' at the equator and 1.6 x 10“^ s“’ at the 
poles. These error ranges are time variable, with uncertainties 
exceeding 2 x 10~^ s“' at the poles in 2005. 

In addition, the quantitative values of Cp (and hence N) de¬ 
pend on an assumption about the equilibration timescale for 
ortho-para-H 2 conversion. The standard assumption is to as¬ 
sume the ‘frozen’ limit, where conversion is extremely slow 
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and the two spin isomers can be treated separately - i.e., Cp in¬ 
cludes individual contributions of ortho- and para-H 2 , in ad¬ 
dition to contributions from He and CH 4 . Latent heat release 
from the conversion between ortho- and para-H 2 is not consid¬ 
ered. Given that disequilibrium is detected throughout much of 
Saturn’s troposphere (Fig. 1 1), a long equilibration timescale 
is a reasonable assumption and will be employed throughout 
this study. However, we cannot completely rule out the alter¬ 
native assumption that equilibration timescales are very short, 
particularly at higher pressures due to catalysis on the surface 
of Saturn’s tropospheric aerosols (Massie and Hunten, 1982). 
Following the review by Irwin (2009), we recalculate the ro¬ 
tational partition function assuming the two extremes: either an 
equilibrium mix or treating the two gases separately. The dif¬ 
ference at equinox can be seen in Fig. 15 - as the equilibrium 
adiabatic lapse rate gfcp is smaller than the ‘frozen’ lapse rate, 
N is systematically smaller by =» 1 x 10“^ s“' throughout the 
upper troposphere. The troposphere is therefore more unstable 
to buoyancy oscillations than in the frozen case, and N becomes 
imaginary for p > 500 mbar (the R-C boundary). This system¬ 
atic offset, in addition to the random uncertainties, limits our 
ability to quantify the buoyancy frequency until some indepen¬ 
dent estimate of Saturn’s static stability is available. 

For the ‘frozen’ assumption, F < Fa throughout the 70-600 
mbar region in Fig. 15, implying that remains positive and 
the upper troposphere is sub-adiabatic and stably stratified, as 
expected. The tropopause has a higher static stability at the 
equator than at the neighbouring belts, and the stability tends 
to increase from equator to pole, following similar patterns to 
those derived from Voyager/IRIS data (Fig. 11 of Conrath et al., 
1998, which had N ^ 9 x 10“^ s“^ at 200 mbar consistent with 
our Fig. 15). Small changes are apparent from 2005 to 2013, 
with the southern hemisphere showing a decreasing N in the 
500-mbar region (from 6 x 10“^ s“* to 4 x 10“^ s“^), indicating 
a tendency towards more unstable conditions that is larger than 
the 1 X 10~^ s“' uncertainty. The stability becomes smaller at 
deeper levels and the equator shows some of the smallest values 
of N near 500 mbar (3 x 10“^ s“^) but this has not changed with 
time. In the equilibrium case, the values of N are lower and the 
R-C boundary is found for p > 500 mbar. The tropical region 
within +30° latitude of the equator shows the most substantial 
changes over time, associated with the warming of the NEB 
and cooling of the SEB, but these remain sub-adiabatic at all 
times. Einally, northern mid-latitudes show a change between 
the 2009 equinox and the present day that is comparable to our 
levels of uncertainty, as the warming due to Saturn’s northern 
storm (Eletcher et al., 2011; Achterberg et al., 2014) decreased 
F and therefore also N. Unfortunately we cannot determine 
whether this stratification change preceded the storm eruption 
(i.e., to relate upper tropospheric stability to the storm trigger¬ 
ing) due to the time-sampling of the CIRS data. These ‘frozen’ 
estimates of N will be used in the calculation of potential vor- 
ticity gradients in the following sections. 

5.2. Tropospheric Winds 

If we assume geostrophic balance, then the shifting merid¬ 
ional temperature gradients {dT jdy, where y is the north-south 


distance) are related to the vertical shear on zonal jets (dujdz) 
via the thermal windshear equation (e.g., Andrews et al., 1987, 
pl20). Read et al. (2009) used Cassini temperature observa¬ 
tions during southern summer (Fletcher et al., 2007) to present 
a snapshot of the vertical windshear, observing a general decay 
of the zonal jets with height that was strongest in the southern 
stratosphere. They called for further measurements of Saturn’s 
wind field throughout the seasonal cycle, so we use the thermal 
windshear equation in the zonal direction: 


du 

(9 In (p) 


Rdl^j^dl^_ d_u 
a dijj dy dz 


( 2 ) 


This equation allows us to integrate the zonal winds u(if/) (where 
ij/ is the latitude) as a function of altitude (z), and therefore to 
explore gradients in the wind fields in the horizontal and ver¬ 
tical directions. / = 2f2sin(i/r) is the Coriolis parameter (the 
vertical component of the planetary rotation vector, where Q is 
the planetary angular velocity); aixjj) is the planetary radius; R 
is the molar gas constant divided by the mean molar weight of 
Saturn’s atmosphere; and H - RT/g is the scale height of the 
atmosphere. 

In this work, the interpolated temperature fields are used to 
calculate the changing dujdz with time (Fig. 16), which is inte¬ 
grated to give u{z,(ji), assuming the wind fields derived from 
continuum-band Cassini imaging by Garcfa-Melendo et al. 
(2011) as a time-invariant boundary condition at 500 mbar. 
Garcfa-Melendo et al. (2011) calculated these velocities using 
the System III rotational reference frame measured by the Voy¬ 
ager mission Seidelmann et al. (2007). Fig. 16a shows a typical 
latitude-altitude cross section of Saturn’s tropospheric winds 
in 2009, and we find that seasonal variations from this struc¬ 
ture are small. As the thermal wind equation involves spa¬ 
tial gradients of the interpolated temperature fields, the MC 
technique was used to investigate the range of uncertainties in 
u{z,\j/), shown in Fig. 16b. Uncertainties are smallest at 500- 
mbar where we specify the cloud-top winds (although these ne¬ 
glect the 5-10 m/s uncertainties on the cloud tracking measure¬ 
ments, Garcfa-Melendo et al., 2011), and grow with altitude 
due to the uncertainties in integration, being of order 10 m/s at 
the tropopause. We caution that the absolute magnitude of the 
zonal winds at each altitude is extremely sensitive to the uncer¬ 
tainty in the location of the cloud-tracked winds, particularly if 
the altitude of the cloud tracers is latitudinally variable, which 
would add further systematic error to our measurements. 

Fig. 16c shows that changes over the ten-year span of these 
observations are small at most latitudes, and within the uncer¬ 
tainty derived from the MC simulation. Variations only appear 
significant in Saturn’s northern hemisphere, particularly for the 
broad retrograde jets in the hemisphere emerging from winter. 
The sharper, narrower prograde jet peaks experience negative 
shear with altitude and appear largely unaffected over time. In 
the southern hemisphere, we see a trend of increasingly posi¬ 
tive (prograde) shear on the broad westward jets from summer 
through to autumn (Fig. 16d), but this has limited effect on 
the zonal winds at 100 mbar (Fig. 16c). In the north, the pic¬ 
ture is more complex. Westward jets at 100 mbar poleward of 
45 °N have become more eastward with time by approximately 
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(a) Brunt Vaisala Frequency [10'^ s'^] 


Planetographic Latitude 

-90 -75-60 -45 -30 -15 0 15 30 45 60 75 90 



(b) Buoyancy Frequency Uncertainty 



Planetographic Latitude 


(c) Alternative Buoyancy Frequency (2009) 



Planetographic Latitude 


Planetographic Latitude 


Figure 15: (a) The Brunt Vdisdld buoyancy frequency (10^^ s“*), a diagnostic of the stratification of Saturn’s upper troposphere, as a function of latitude and pressure 
for Saturn’s northern spring equinox ±4 years, calculated from the temperature and gravity field assuming ‘frozen’ equilibrium of ortho-para-H 2 as described in 
the main text. Spurious oscillations at high north polar latitudes in 2005 are omitted, as these were not directly constrained by CIRS observations, (b) Uncertainty 
range on N (10“^ s“^) in 2009 as a representative example of the time series, calculated via the MC simulation, (c) Alternative estimate of the buoyancy frequency 
assuming instantaneous equilibration between the two spin isomers of H 2 leading to a generally more unstable troposphere. 


19 










































10 m/s, corresponding to an increasingly positive (prograde) 
vertical shear over time. Although this change is equivalent to 
our uncertainty, the trend with time is continuous rather than 
random. The exception is the westward jet at 39°N (associ¬ 
ated with the outbreak of the 2010 disturbance), which has be¬ 
come more westward with time due to the increasingly retro¬ 
grade (negative) vertical shear. This is consistent with the post¬ 
storm evolution of the thermal field noted by Achterberg et al. 
(2014) and the cloud-tracked wind velocities from Sayanagi 
et al. (2013). 

The classical thermal wind equation relies on an evaluation 
of the Coriolis parameter, and therefore becomes unreliable at 
low latitudes, as indicated by the growth of errors in Fig. 16. 
Nevertheless, the flanks of the equatorial jet show some inter¬ 
esting variability, with the vertical shear at 10°N becoming in¬ 
creasingly negative (retrograde) over time. The lack of similar 
changes at 10°S suggests a process acting in the northern hemi¬ 
sphere that is absent from the south. Interestingly, this process 
is making the zonal windshear profile in Fig. 16d appear more 
symmetric about the equator with time, and this latitude (10°N) 
was identified by Li et al. (2013) as the location of the largest 
difference between tropospheric winds derived from Voyager 
and Cassini data. This difference between the hemispheres 
could be caused by a cross-equatorial flow such as those hy¬ 
pothesised by Friedson and Moses (2012) for higher altitudes. 
Alternatively, interactions between the zonal-mean flow and 
waves emitted from the northern storm and ‘beacon’ (Fletcher 
et al., 2012) might also serve to modify the flow of the equa¬ 
torial jet on its northern flank. In summary, while the southern 
summer/autumn jet system appears to evolve very slowly, the 
changes to the northern jets are more substantial, possibly as a 
result of (i) the downward transport at high latitudes suggested 
by the para-H 2 distribution, or (ii) the seasonal heating of the 
upper troposphere and associated hazes; or (iii) the disruption 
caused by Saturn’s northern storm in 2010-11. 

5.3. Tropospheric Potential Vorticity 

Having derived Saturn’s tropospheric wind held u(z, ij/), we 
now explore the implications of gradients in these zonal Hows 
on dynamic activity. Our goal is an understanding of the possi¬ 
bilities for wave propagation in the upper troposphere, but we 
start by assessing quasi-geostropic potential vorticity (QGPV, 
qc) gradients as a function of latitude, altitude and time. QGPV 
gradients have previously been studied for Jupiter and Saturn 
(Read et al., 2006; Read et al., 2009), showing them to be a 
powerful diagnostic tool for understanding dynamic Hows in a 
stably-stratiHed atmosphere. We assume that quasi-geostrophy 
is valid given the small Rossby numbers associated with Sat¬ 
urn’s atmospheric flows (i.e., Rossby numbers smaller than 
one). The meridional QGPV gradient (sometimes called ‘ef¬ 
fective beta,’ Pe - dqddy, e.g., pl27 of Andrews et al., 1987) 
is assembled from the sum of three terms: (1) p - dfidy, 
the northward gradient of planetary vorticity (i.e., the Corio¬ 
lis parameter); (2) yS,. = -d^ujdy^, the meridional curvature 
of the zonal wind Held; and (3) yS^, the vertical curvature of 
the wind field. This Pe is only strictly valid when the merid¬ 
ional length scales are smaller than the planetary radius and 


longitudinal temperature contrasts are small. Both approxima¬ 
tions hold for the scales of Saturn’s belt/zone motions consid¬ 
ered here. Following the definitions of Andrews et al. (1987); 
Sanchez-Lavega et al. (201 1): 


Pe - P+ Py -H Pz 

d^u I d I f du\ 
^ dy^ p dz \ dz j 


(3) 

(4) 


Each variable has been defined previously, with the exception of 
the atmospheric density p - pIRT from the ideal gas equation. 
The full two-dimensional (latitude, pressure) variation of each 
term (gravity, heat capacity, temperature, scale height, density, 
zonal wind, buoyancy frequency) is accounted for in the cal¬ 
culation. Py therefore uses the zonal wind field calculated at 
all altitudes via the thermal windshear relation (Section 5.2); Pz 
uses the vertical windshear and buoyancy frequency (Section 
5.1). 

The effective p and its three components are shown in Fig. 
17(b) as a function of latitude for 2009 at the 330-mbar level 
(i.e., averaging the reconstructed temperature fields and winds 
for the year surrounding Saturn’s equinox). Values near the 
equator are omitted, as the small Coriolis parameter limits the 
reliability of the extrapolated wind field at low latitudes. The 
Py and Pz terms introduce a strong variability of the QGPV gra¬ 
dient with latitude, although the two are largely anti-correlated 
with one another. Unfortunately, the necessity for taking spatial 
gradients of the temperatures, winds and para-H 2 fields serves 
to magnify uncertainties, so we use the MC simulation (100 
realisations of Gaussian perturbations to the 3000-1- T{p) and 
fpip) profiles) to quantify the range of errors on yS,,, ySj, and Pg in 
Fig. 17c. As expected, the uncertainty on Pg (shown for 2009 
in Fig. 17d) is dominated by the vertical curvature of the wind 
field Pz because of the trend of increasing uncertainty on the 
zonal wind with height. In comparison, the meridional cur¬ 
vature of the winds Py provides a smaller contribution to the 
uncertainty, given that the meridional gradients of uncertainty 
in the winds in Fig. 16b are rather small. Fig. 17d also shows 
where the uncertainty in Pg becomes so large (> 20 x 10“'^ 
m 's ') that it dwarfs the meridional changes in Pg, i.e. for 
p > 600 mbar in the mid-latitudes and for p > 300 mbar at 
high latitudes. QGPV gradients are more uncertain in the north¬ 
ern hemisphere due to the sparser time series and cooler winter 
temperatures (i.e., larger retrieval uncertainties). 

Despite the large uncertainties. Fig. 17b-c indicate that Pg 
shows reversals in sign as a function of latitude that are larger 
than the errors, particularly equatorward of +60°. QGPV gra¬ 
dients at high latitudes are too uncertain to draw conclusions. 
The strongest gradients are associated with the peaks of the pro¬ 
grade jets, as previously discovered by Read et al. (2009). The 
QGPV gradient reversals result in a ‘staircase pattern’ in the po¬ 
tential vorticity, with steep steps (i.e., high gradients) associated 
with the prograde jet peaks, and a tendency to flatten in between 
(i.e., homogenising the QGPV) in association with minima in u. 
Read et al. (2009) suggested that the presence of strong positive 
QGPV gradients indicates that the jets are dynamical bound¬ 
aries at the edges of regions of intense latitudinally-localised 
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(a) 2009 Thermal Wind [m/s] 
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(c) Zonal Wind Change at 100 mbar 



(d) Vertical Wind Shear at 100 mbar 
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Figure 16: Zonal mean winds predicted by integrating the thermal windshear equation using the cloud-top winds of Garcia-Melendo et al. (2011). Winds are 
relatively stable, with panel (a) showing the annual-average zonal winds in 2009 (equinox) with a contour spacing of 30 m/s - regions within ±10° of the equator 
are omitted because the thermal wind balance becomes invalid at low latitudes. Panel (b) shows the range of uncertainty on the zonal winds in 2009 from the 
MC simulation with contours spaced every 2 m/s. Panel (c) shows the small changes in the zonal wind field at 100 mbar (solid lines for different values of L.y) 
compared to the uncertainty envelope (dotted line). Panel (d) shows the zonal windshear per scale height at 100 mbar for the same six Ls values, also compared to 
the uncertainty envelope (dotted line). 
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Figure 17: Evolution of Saturn’s tropospheric quasi-geostrophic potential vor- 
ticity (QGPV) gradients as a function of time and latitude. Panel (a) shows 
the cloud-top zonal winds from Garcia-Melendo et al. (201 1), comparing both 
continuum-band (black) and methane-band (red) measurements of the wind 
speeds which sample deep winds (k 500 mbar) and upper tropospheric winds, 
respectively. Panel (b) shows the contributions to the ‘effective beta’ (jSe) 
QGPV gradient at 330 mbar at Saturn’s equinox in 2009, and panel (c) shows 
the uncertainties in these values from the MC simulation. Panel (d) shows the 
variation of uncertainty with latitude and pressure, in units of 10“*^ m“*s“* 
- contours higher than 20 X 10“'^ m“*s“* are not shown. Panel (e) shows the 
evolution of fj,, as a function of year (solid lines with colours shown by the inset 
key), which causes changes in the QGPV staircase structure (Read et al., 2009) 
over time. Dotted lines show the time-dependent uncertainty in for 2005, 
2009 and 2013 from the MC simulation. 


activity and mixing of vorticity. The mixing serves to remove 
the initial PV gradients (flattening), moving the sharpest gradi¬ 
ents to the edges of the band to be associated with a local jump 
in QGPV across the prograde jets. 

A sign change for the PV gradient with latitude is a necessary 
(but not sufficient) condition for instability with respect to the 
Charney-Stern condition (Charney and Stern, 1962). Read et al. 
(2009) further point out that Saturn’s most interesting dynamic 
phenomena seem to occur near to where the QGPV gradient 
changes sign - the polar hexagon, ribbon wave and ‘storm alley,’ 
and that these are locations where the mean zonal flow appears 
to violate the neutral-stability criterion relative to Arnol’d’s sec¬ 
ond theorem. A similar sign reversal is evident in the gradi¬ 
ent of absolute vorticity (fi - d^ujdy^) (Sanchez-Lavega et al., 
2000), but note that a large proportion of the QGPV gradient 
reversals are driven by the vertical curvature term (ySj;, some¬ 
times referred to as the ‘stretching’ component to distinguish it 
from the gradient of absolute vorticity), and hence the retrieved 
temperature field from CIRS is essential for this analysis. How¬ 
ever, a more thorough analysis of the relation between the mean 
zonal flow and stability violation will likely require numerical 
modelling of the fluid flows. 

The seasonal trends in temperature identified in this work 
permit an assessment of the variability of the ‘staircase struc¬ 
ture’ in QGPV, as shown in Fig. 17(e). The QGPV gradi¬ 
ents have evolved over the ten years since the analysis of Read 
et al. (2009), but we caution that the uncertainty on fie is time- 
dependent, being larger at the start and end of the time series 
than in 2009. Annual-average uncertainties are depicted as 
dotted error envelopes in Fig. 17e, and show that fie changes 
with time are dwarfed by the uncertainty. Despite this quan¬ 
titative uncertainty, qualitative trends are visible: at southern 
latitudes, we observe a weakening of the QGPV gradient as¬ 
sociated with the 60.9° S prograde jet, whereas those associ¬ 
ated with jets at 73.9°S and 48.2°S have remained largely un¬ 
altered. In between the 60.9°S and 73.9°S jets, the QGPV gra¬ 
dient has shifted from being very negative to being flatter, and 
more ‘step-like.’ Tropical latitudes between 30°S and 30°N, 
dominated by Saturn’s broad equatorial flow, have remained 
largely unaltered over time. The picture in the northern hemi¬ 
sphere is extremely confusing: uncertainties poleward of 60°N 
are so large that the changes are not deemed trustworthy there. 
However, the QGPV gradient associated with the prograde jet 
at 47.8°N has altered substantially throughout the mission, out¬ 
side of the uncertainty envelope. We note that this jet marked 
the edge of the ‘warm anomaly’ identihed at depth in the north¬ 
ern winter hemisphere by Fletcher et al. (2007) (Fig. 6a) that 
disappeared as the mission progressed, and strong QGPV gra¬ 
dients were already noted here by Read et al. (2009). We spec¬ 
ulate that the changing sign of the QGPV gradient could have 
affected the stability properties of the atmosphere in this region 
prior to the storm eruption in 2010 in the nearby retrograde jet 
at 39°N. Furthermore, the storm eruption near 39°N could have 
provided intense vorticity mixing that pushed the sharpest gra¬ 
dients of QGPV northward to the prograde jet at 47.8°N. Firm 
conclusions about the stability properties in Saturn’s northern 
mid-latitudes would require a more focussed study of QGPV 
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gradients surrounding the storm epoch. 

In summary, the meridional gradient of QGPV has been as¬ 
sessed as a function of time, latitude and pressure, and was 
found to be variable over the ten year span of Cassini’s mea¬ 
surements despite the large uncertainties involved in this cal¬ 
culation. The implications for these gradient changes, both in 
terms of the vigour of atmospheric mixing and the violation of 
stability criteria, remain poorly understood, but QGPV derived 
from temperatures and wind data may prove to be an invaluable 
tool for studying giant planet meteorology. Finally, the chang¬ 
ing QGPV gradients could potentially act as a baroclinic source 
of planetary waves (Achterberg and Flasar, 1996), whose prop¬ 
agation we discuss in the next section. 


5.4. Planetary Wave Activity 

A final application of the CIRS-derived seasonal tempera¬ 
ture profiles and associated QGPV gradients is a study of po¬ 
tential wave activity in Saturn’s troposphere. Achterberg and 
Flasar (1996) previously used Voyager/IRIS temperature re¬ 
trievals to show how the meridional and vertical curvature of 
the wind fields served to confine the propagation of planetary 
waves (or Rossby waves) to particular latitude and altitude do¬ 
mains. Rossby waves owe their existence to the change in the 
Coriolis force with latitude acting as a restoring force (i.e., /?«), 
but the ability to propagate in the vertical depends on the in¬ 
dex of refraction of the background medium (v^) (Chamey and 
Drazin, 1961), which can be estimated from the quantities de¬ 
rived above. The dispersion relation for a 3D Rossby wave in 
a baroclinic atmosphere (assuming a beta plane) allows us to 
estimate the zonal phase speed of the wave (c^) relative to the 
mean zonal wind (m) as follows: (Sanchez-Lavega et al., 201 1): 

_ _ _ 

“ ^"“^2 + ^2 + (/ 2 / A 2 )( m 2 + n 2 ) 

where k, I and m are the zonal, meridional and vertical wave 
numbers of the wave, respectively, and n^ = l/( 4 // 2 ) following 
Sanchez-Lavega et al. (201 1). However, this omits higher-order 
vertical derivatives of the buoyancy frequency (Achterberg and 
Flasar, 1996), which are important for Saturn as N (i.e., the 
static stability) is increasing from near zero in the convective 
troposphere to large stratospheric values in Fig. 15. We in¬ 
clude these additional terms in n^ following Karoly and Hoskins 
(1982): 

1 ...... 
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where the last equality re-expresses the derivatives in terms of 
pressure. Following Achterberg and Flasar (1996), we rewrite 
the dispersion relation into a form where parameters related to 
the background medium (m, N'^, H, f and Pe) are separated from 
the characteristics of the wave: 
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where v is the index of refraction that depends only on the back¬ 
ground atmosphere and the phase speed of the wave. The final 
equation uses the pressure-coordinate version of n^ to provide a 
useable expression for v2 based on the parameters of the back¬ 
ground atmosphere. If y2 > 0 then the refractive index is real, 
permitting wave solutions that can propagate. Alternatively, 
v2 < 0 implies an imaginary refractive index and evanescent 
waves that can decay exponentially. The boundary at v2 = 0 
can serve as a trapping barrier, keeping the wave within the 
region of real refractivity via reflections at this boundary (al¬ 
though there will be leakage through the region of imaginary 
refractive index). Note that when u - Cx = 0, v will tend to¬ 
wards infinity - these are known as critical surfaces that could 
also act to confine the waves, causing wave absorption, wave 
breaking and interactions with the mean zonal flow (Achter¬ 
berg and Flasar, 1996), although some transmission of energy 
through this critical surface remains possible. Here the wave 
solutions become non-linear and unsteady, and can lead to the 
creation of vortices. In summary, an estimate of v2 based on 
the properties of the background medium permits a crude ex¬ 
ploration of the potential for Rossby wave activity in Saturn’s 
troposphere. 

However, this derivation of the index of refraction is a sim¬ 
plification, as it assumes separability of the wave solutions 
in the zonal, meridional and vertical dimensions. Here we 
use the Wentzel-Kramer-Brillouin-Jeffries (WKBJ) approxima¬ 
tion (Section 4.5 of Andrews et al., 1987), whereby the wave¬ 
lengths are assumed to be much smaller than the length scales 
over which the background atmospheric medium varies. This 
will not necessarily be the case in the real atmosphere (par¬ 
ticularly for small wavenumbers), and wave solutions will not 
be sinusoidal as they will vary spatially in the meridional (/) 
and vertical (m) directions. We have implicitly assumed that 
the wave solution is locally sinusoidal in Eq. 5. A more com¬ 
prehensive use of the WKBJ approximation would require the 
launching of rays in multiple directions to properly compute the 
meridional and vertical wavenumbers from numerical raytrac¬ 
ing. Nevertheless, the equations above can still provide useful 
basic insights into the propagation of planetary waves in the 
middle atmosphere (Andrews et al., 1987). 

Using the two-dimensional fields of fie, N^, H, u and /, we 
calculate v2 as a function of time, latitude and pressure in the 
left column of Fig. 18. For simplicity we assume that Cx = 0 
(i.e., phase speeds stationary with respect to the background 
flow), which implies that critical surfaces only exist for u ^ 0 
m/s in the 2D zonal wind held. This zero-zonal-velocity surface 
(i.e., the critical surface where v —> oo for stationary waves) 
is included as blue lines in Fig. 18. Regions where v2 > 0 
(i.e., a real refractive index permitting wave propagation) are 
shown in grey, compared to those regions where v2 < 0 (i.e., an 
imaginary refractive index implying evanescent regions), but 
we warn the reader that the vertical extent of the positive-index 
regions will be sensitive to the assumption of Cx = 0. Both the 
critical surfaces and the v2 = 0 boundary appear to change with 
time in Fig. 18. 

The calculated values of v2 rely on spatial derivatives of the 
temperature, wind and composition held (via the calculation of 
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N^), and are therefore subject to large uncertainties that grow 
extremely rapidly with depth. The MC simulation permits a 
crude estimate of the range of expected values based on the 
temperature and para-H 2 uncertainties, and annual-average un¬ 
certainties are shown in the right-hand column of Fig. 18. Re¬ 
gions of the figure where uncertainties exceed 1x10“'^ nT^ (the 
contour spacing on the figure) have been shaded to indicate 
low conhdence in the results, and these coincide with regions 
where the zonal windheld u approaches zero (i.e., the critical 
surfaces). As a result, our knowledge of the atmospheric refrac- 
tivity of Rossby waves is extremely sensitive to uncertainties in 
Saturn’s winds. The most trustworthy regions are in the trop¬ 
ics, within approximately +35° of the equator and for p < 400 
mbar. 

The scale of this uncertainty on is rarely reported in the 
literature, and precludes detailed quantitative analysis. How¬ 
ever, some qualitative trends are apparent that will be worthy 
of future analysis. In the specific case of stationary waves, real 
refractive indices appear to be conhned to p < 300 mbar by 
the static stability. These > 0 regions are found surrounding 
the equator in the 100-200 mbar region with latitudinal bound¬ 
aries that are relatively stable, but vertical extents than change 
with time. For example, the positive-index region near 10-15°N 
extends to lower pressures as time progresses, and the positive- 
index region near 15 - 25°S shrinks in extent over time. Ad¬ 
ditional regions of positive occur at mid-latitudes in both 
hemispheres, but here the uncertainties on have increased 
substantially. The vertical extension of these positive regions 
appears to change from year to year, but trends could be spuri¬ 
ous given the level of uncertainty. We compare our results with 
Fig. 8 of Achterberg and Flasar (1996), who used ray tracing 
calculations to understand the propagation of quasi-stationary 
waves based on inversions of Voyager/IRIS spectra. They also 
found regions of real refractivity, one equatorward of 25°N that 
moved closer to the equator with increasing depth (a similar be¬ 
haviour is seen in Fig. 18); and a branch of positive refractivity 
between 30 - 40°N that extended down into the upper tropo¬ 
sphere, just as we are finding with the Cassini/CIRS data. They 
concluded that the propagation of quasi-stationary waves was 
limited to these latitudes by the curvature of the wind field (i.e., 
jSj), and limited to p < 200 mbar by the static stability (i.e., /3^), 
in broad agreement with our hndings from Cassini/CIRS. 

Adding further conhdence to the distributions in Fig. 18, 
the southern mid-latitude region (35 -45°S) and northern equa¬ 
torial region (0-30°N) have both been identihed as locations of 
strong stratospheric wave activity during the Cassini epoch (Or¬ 
ton et ah, 2014), implying a connection to the refractivity pro- 
hles of the upper troposphere. Furthermore, a near-stationary 
thermal wave (wavenumber 2) was identihed by Voyager IRIS 
at 130 mbar at northern mid-latitudes during spring (Achterberg 
and Flasar, 1996). In the present analysis, we show that the ver¬ 
tical domain of real refractivity varies because of seasonal shifts 
in the static stability (e.g.. Section 5.1), whereas the curvature 
of the wind held is less variable so that the broad meridional 
structure of the refractive index remains the same over time. 
We speculate that slow seasonal changes to the thermal struc¬ 
ture, and hence the propagation of planetary waves, inhuences 


the variable storm activity with time; for example, the southern 
‘storm alley’ (38.5-46°S planetographic, Vasavada et al., 2006) 
is found to be most active during southern summer/autumn; 
whereas the eruption of Saturn’s northern storm and its inhu- 
ence on the stratosphere (Fischer et al., 2011; Fletcher et al., 
2011) occurred during northern spring. As wave activity ap¬ 
pears to be commonplace in Saturn’s stratosphere (e.g. Orton 
and Yanamandra-Fisher, 2005), and this is where waves are 
likely to break and deposit their energy (Achterberg and Flasar, 
1996), an extension of this study to lower pressures is required, 
in addition to proper diagnostics using k, I, m and Cx for ob¬ 
served planetary waves on Saturn (Orton et al., 2014). 

6. Conclusions 

Far-infrared Cassini/CIRS spectra of Saturn’s hydrogen and 
helium continuum have been used to track the seasonal evo¬ 
lution of the upper troposphere over a third of a Saturnian 
year (2004-2014). Inversions of CIRS spectra allowed us to 
determine three-dimensional prohles (time, latitude, pressure) 
of temperatures, para-H 2 , zonal winds, stratihcation, potential 
vorticity gradients and the refractive index for wave propaga¬ 
tion. This builds upon earlier work that captured a snapshot 
of tropospheric conditions during southern summer (Fletcher 
et al., 2007), and allows a quantitative comparison between 
Voyager (1980) and Cassini (2010) observations at the same 
point in Saturn’s seasonal cycle. Our conclusions are as fol¬ 
lows: 

1. Thermal Asymmetries: Wintertime asymmetries in 

tropospheric temperatures have largely vanished by late 
spring, so that temperatures are approximately symmetric 
about the equator. The timing of the northern tempera¬ 
ture minima and southern temperature maxima varied ac¬ 
cording to latitude and altitude (occurring later at higher 
pressures), but summertime maxima occurred sooner after 
solstice than the wintertime minima, some of which oc¬ 
curred after the spring equinox (2009). These trends are 
qualitatively consistent with the expectations of radiative- 
convective simulations, but quantitative comparisons are 
hampered by the uncertain heating effects within the up¬ 
per tropospheric hazes. Saturn’s upper troposphere has re¬ 
mained sub-adiabatic and stably-stratihed (positive buoy¬ 
ancy frequencies) throughout the period of observations 
(assuming slow interconversion between ortho- and para- 
H 2 ), confirming that CIRS senses altitudes above the con¬ 
vective troposphere. Contrasts due to Saturn’s belt/zone 
structure are superimposed onto the seasonal trends, and 
maxima in dTfdy remain co-located with Saturn’s cloud- 
top zonal jets. 

2. Aerosol Heating: The thermal perturbation in the 100- 
250 mbar region due to localised heating in Saturn’s haze 
layer has evolved over time. The magnitude of the T{p) 
curvature fell over the southern hemisphere (vanishing 
poleward of 45°S by 2014) and grew over the northern 
hemisphere (starting with the cool zone near 45°N after 
2008) as spring progressed. By 2014, the NEB and SEB 
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Figure 18: Left column: Three examples of annual-average zonal mean contours of the index of refraction (v^ in 10“^^ m““) for planetary waves in Saturn’s upper 
troposphere, calculated for the specific case of stationary waves (Cx = 0). Contours greater than ±1 x 10“*^ m~^ are shown every 1 x 10“^^ m~^, with dense contours 
suggesting rapidly increasing v close to critical surfaces. Contours less than ±1 x 10“^^ m“^ are shown every 0.25 x 10“^^ m~^. Regions of real refractive index 
(v" > 0) where planetary waves may propagate are shaded grey and feature solid-line contours. Regions of imaginary refractive index (v^ < 0), where waves are 
evanescent, are white and feature dotted contours. Critical surfaces are represented by the blue curves where zonal-mean winds are zero (i.e., where v —> oo). The 
equatorial region is omitted due to uncertainties in the wind field. Right column: Uncertainties in (10“^^ m“^) calculated via the MC simulation for the same three 
years. Uncertainties smaller than 1 x 10“^^ m“^ (i.e., the contour spacing in the left column) are shown with contours spaced by 0.02 x 10“^^ m“^. Uncertainties 
larger than 1 x 10“^^ m“^ are shaded grey and have contours spaced every 5 x 10“^^ 
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exhibited thermal perturbations of the same relative size. 
Unfortunately, CIRS data alone cannot assess whether the 
shifting magnitude of the perturbation is due simply to 
enhanced sunlight on a pre-existing and static haze (e.g. 
Friedson and Moses, 2012), or whether the dominant ef¬ 
fect is due to the changing colour and opacity of the hazes 
observed by Cassini (e.g., Edgington et al., 2012). 

3. Dynamic Heating: Anomalously warm regions in the 
deep troposphere (p > 500 mbar) evolve over non- 
seasonal timescales. These include a warm region between 
50-70°N first identified by Fletcher et al. (2007), the heat¬ 
ing due to slow subsidence in the aftermath of Saturn’s 
springtime storm at 25 - 45°N (Achterberg et al., 2014), 
and a similar warm patch identified in Voyager retrievals 
between 25 - 45°N (this work and Conrath and Pirraglia, 
1983). These warm regions are at the high-pressure limit 
of CIRS vertical sensitivity, but could be a dynamical re¬ 
sponse to the changing stratification of the overlying at¬ 
mosphere. 

4. Para-H 2 : Changes in the para-H 2 fraction with time are 
subtle and influenced by dynamics, fp showed signs of a 
general increase in the northern hemisphere (0.02-0.03 at 
200-500 mbar) when the strong effects of the low-/^ air 
lofted by the 2010 storm are discounted. Northern lati¬ 
tudes poleward of 60°N show the largest increases near 
the tropopause, whereas the north polar fp at 500 mbar 
was largely unchanged. This has the effect of changing the 
symmetric fp distribution of 2005 into a more asymmetric 
distribution in 2014 (with higher fp in the north), sugges¬ 
tive of general subsidence in the springtime hemisphere. 
This brings the CIRS 2010 results into line with the Voy¬ 
ager 1980 results, although the latter showed a stronger 
asymmetry suggesting that para-H 2 variability (and there¬ 
fore the upper tropospheric subsidence) does not proceed 
at the same rate every Saturnian spring. 

5. Para-H2 Disequilibrium: Disequilibrium in para- 

H 2 reflects primarily the changing temperature structure 
(and therefore distribution of fp^eqm) rather than actual 
changes in fp induced by chemical conversion or transport. 
The reversal of the para-H 2 disequilibrium asymmetry, 
where strong southern super-equilibrium conditions and 
strong northern sub-equilibrium conditions have both been 
eroded away, is mostly due to the temperature changes and 
not to the para-H 2 changes described in the previous con¬ 
clusion. Neither large-scale interhemispheric transport, 
nor enhanced catalysis on seasonal hazes, are required to 
explain this shifting asymmetry. Nevertheless, dynamics 
can have a strong influence, with (i) up welling creating 
sub-equilibrium conditions at the equator that have de¬ 
clined slightly over time; (ii) upwelling of low-/p air by the 
northern storm (which is still in the process of returning 
the the ‘normal’ conditions at the time of writing); and (iii) 
increasing subsidence in the spring hemisphere rapidly re¬ 
ducing the sub-equilibrium conditions found there during 
winter. 

6. Tropical Oscillations: Saturn’s upper troposphere in the 
tropics is influenced by the wave train linked to the strato¬ 


spheric ‘semi-annual’ oscillation (which should be cor¬ 
rectly referred to as a quasi-periodic oscillation). Equa¬ 
torial tropopause temperatures have fallen 3 + 1 K over ten 
years as a cool anomaly propagated downwards from the 
stratosphere, and the aerosol heating perturbation appears 
to have been pushed deeper from 80 to 150 mbar as time 
progressed. Temperature, aerosol heating and para-H 2 
contrasts in the tropics appear markedly different between 
the Cassini era and the Voyager era (although the general 
hemispheric asymmetries are well reproduced), suggest¬ 
ing that the influence of the stratospheric oscillation is not 
in phase with the seasonal cycle at the tropopause (i.e., it is 
not ‘semi-annual’ at these depths). Indeed, Sinclair et al. 
(2014) also suggested that Voyager and Cassini saw the 
oscillation at a markedly different phase despite being in 
the same season. 

7. Zonal Wind Variability: Variable meridional temper¬ 
ature gradients imply that the vertical shear on Saturn’s 
zonal jets changes with time, by values typically smaller 
than 10 m/s per scale height at the tropopause. The sharper 
eastward jets experience negative shear with altitude and 
are largely unaffected with time; the broader westward 
jets in the southern hemisphere have experienced increas¬ 
ingly prograde shear with as summer gave way to autumn; 
and the general picture in the northern hemisphere was 
complicated by the influence of the springtime storm (the 
westward jet at 39°N has strengthened due to the increas¬ 
ingly retrograde (negative) vertical shear associated with 
the storm heating). Intriguingly, the northern flank of the 
equatorial jet has experienced the largest change in wind- 
shear, whereas the southern flank remained largely unal¬ 
tered, and it is notable that Li et al. (2013) observed that 
largest contrasts between Voyager and Cassini-era winds 
at this latitude (10°N). 

8. Potential Vorticity Variability: The meridional quasi- 
geostrophic potential vorticity (QGPV) gradient (‘effec¬ 
tive beta’), and hence the ‘staircase’ of potential vorticity 
identified by Read et al. (2009), has varied over the course 
of the mission. Tropical latitudes between 30°S and 30‘’N 
have remained largely unaltered over time, and quantita¬ 
tive conclusions at high latitudes are precluded by the large 
uncertainties there. Gradients associated with mid-latitude 
jets have varied, particularly in the northern hemisphere 
associated with both the storm eruption and the disappear¬ 
ance of a ‘warm anomaly’ at high pressures poleward of 
45°N. These changes in Pe have implications for the vigour 
of atmospheric mixing, and the magnitude of the gradient 
reversals (i.e., sign changes for pf) could imply different 
degrees of violation of atmospheric stability criteria. 

9. Rossby Wave Refractivity: Shifting potential vorticity 
gradients influence the locations of real refractive index 
for the propagation of planetary waves, and this may mod¬ 
ify the confinement regions for wave activity. The vertical 
domain of real refractivity is expected to vary because of 
seasonal shifts in the temperature profile and static stabil¬ 
ity, whereas the curvature of the wind held is less variable 
so that the broad meridional structure of the refractive in- 


26 


dex remains the same over time. We find that stationary 
Rossby waves would be able to propagate in broad regions 
near the equator (20°N to 20°S) and 35 -45° in both hemi¬ 
spheres for p < 300 mbar (limited by the static stability), 
consistent with the sparse observations of wave activity 
in Saturn’s upper troposphere and stratosphere. However, 
the uncertainties on the refractive index grow extremely 
large with depth and latitude, so the confinement regions 
for planetary wave activity remain to be tested by further 
observations. 

We aim to continue Cassini’s tropospheric characterisation 
beyond 2014, through to Saturn’s northern summer solstice 
(Lj = 90°) in 2017, and beyond from ground- and space-based 
observatories. The implications of the changing wind fields and 
stratification on the atmospheric stability and the potential for 
wave propagation remain to be thoroughly explored, and we 
hope that improved radiative-convective models (fully account¬ 
ing for aerosol heating) will ultimately allow us to use the devi¬ 
ations from thermal equilibrium to diagnose vertical velocities 
in the troposphere. The analysis presented in this paper shows 
the utility of thermal infrared temperature and para-H 2 sound¬ 
ing in revealing the circulation and seasonal variability in a gi¬ 
ant planet atmosphere, and the importance of a multi-decadal 
dataset in understanding the slow variability on these distant 
worlds. 
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